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PREFACE 


The  Electro- Optical  Systems  Atmospheric  Effects  Library  (EOSAEL)  was 
developed  to  assist  the  characterization  of  the  battlefield  environment.  One 
component  of  that  characterization  focuses  on  the  quantification  of  the  scattering 
properties  of  various  natural  boundary-layer  aerosols,  battlefield  dusts,  and 
inventory  smoke  munitions.  This  report  documents  the  improved  aerosol 
scattering  phase  property  data  incorporated  in  the  1992  version  of  EOSAEL. 

The  aerosol  phase  function  database  (PFNDAT)  is  a  continuation  of  previous 
phase  function  databases  released  with  the  1980,  1982,  and  1987  versions  of 
EOSAEL  (c.f.,  Duncan  (ed.),  1980;  Shirkey  et  ah,  1987). 
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ABSTRACT 


The  Electro-Optical  Systems  Atmospheric  Effects  Library  (EOSAEL)  Phase  FunctioN 
DATabase  (PFNDAT)  consists  of  a  series  of  phase  functions  and  extinction  and  scattering 
coefficient  data  for  30  naturally  occurring  and  8  manmade  aerosols  associated  with  the  near 
surface  atmosphere.  These  phase  functions  are  useful  in  characterizing  the  near  surface 
atmosphere  for  propagation  and  scattering  studies  where  typical  scattering  species  are 
required.  Models  using  this  database  include  several  EOSAEL  modules.  The  naturally 
occurring  aerosols  consist  of  the  maritime,  urban,  and  rural  aerosol  size  distributions  at 
eight  relative  humidities  each,  two  fog  distributions,  three  rain  distributions,  and  one  snow 
distribution.  The  manmade  aerosols  consist  of  three  dust  types  and  five  smoke  types. 
The  dusts  include  a  high-explosive  dust  distribution  and  light  and  heavy  loading  dust 
distributions.  The  smokes  consist  of  white  phosphorus  results  for  three  different  relative 
humidities;  hexachloroethane;  and  fog  oil  smoke.  The  database  includes  information  at 
a  variety  of  wavelengths  for  each  scattering  species  (dependent  on  availability  of  index 
of  refraction  data).  Wavelengths  range  from  0.35  to  40.0  fim.  This  report  contains  brief 
descriptions  of  the  aerosol  size  distribution  characteristics,  aerosol  index  of  refraction  data 
used  to  generate  the  phase  functions,  information  on  the  contents  of  the  PFNDAT  database, 
and  graphs  of  the  phase  functions.  This  version  of  PFNDAT  improves  on  the  original 
database  by  increasing  the  resolution  of  visible  band  phase  function  results.  Previous 
versions  only  included  a  0.55  fj,m  result.  Improvements  in  the  AGAUS  code  used  to  generate 
the  Mie  scattering  results  used  has  also  resulted  in  more  accurate  phase  functions. 
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EXECUTIVE  SUMMARY 

Introduction 

The  Aerosol  Phase  Function  Database  (PFNDAT)  is  the  primary  repository  of 
aerosol  scattering  property  data  provided  within  the  Electro-Optical  Systems 
Atmospheric  Effects  Library  (EOSAEL).  EOSAEL  is  designed  to  provide  detailed 
atmospheric  effects  models  of  the  battlefield  environment.  Often  these 
environmental  effects  include  the  influences  of  extinction,  absorption,  and 
scattering  due  to  aerosols.  The  PFNDAT  database  provides  the  necessary 
information  to  characterize  a  wide  variety  of  battle- induced  and  natural  aerosols. 
Included  in  this  collection  are  data  sets  describing  the  three  main  classes  of 
hazes  (rural,  maritime,  and  urban)  at  varying  relative  humidity  levels.  Two 
types  of  fog  (radiation  and  advection)  are  provided.  Rain  properties  are 
provided  at  three  precipitation  rates  (drizzle  (1  mm/h),  moderate  (5  mm/h),  and 
heavy  (10  mm/h)).  One  class  of  snow  is  considered,  along  with  several  battlefield 
induced  contaminants:  Dusts  are  treated  under  three  categories:  light,  heavy 
advection,  and  high  explosive  munition  caused.  Inventory  smokes  include  white 
phosphorus  at  three  relative  humidities,  hexachloroethane,  and  fog  oil. 

Purpose 

Providing  scattering  information  in  a  database  format  allows  aerosol  scattering 
species  to  be  used  within  larger  radiative  transfer  models  without  large  runtime 
calculational  costs.  The  Mie  scattering  codes  used  to  assess  the  scattering 
properties  of  these  aerosols  require  long  processing  times.  When  preprocessed, 
the  resulting  propagation  characteristics  can  readily  be  used  within  the  radiative 
transfer  codes  under  a  variety  of  conditions  without  again  incurring  these  costs. 
Several  EOSAEL  modules  require  input  from  PFNDAT  to  perform  radiative 
transfer  calculations.  Also,  aerosols  are  often  mixed,  requiring  a  weighted 
addition  of  multiple  scattering  species. 

Overview 

Previous  versions  of  PFNDAT  have  focused  on  characterizing  scattering 
properties  within  infrared  (IR)  bands  associated  with  midband  and  far  IR 
sensor  systems.  The  current  calculations  extend  the  previous  results  to  a  higher 
resolution  in  the  visible  band,  and  extend  the  IR  calculations  to  40  /xm.  This 
allows  PFNDAT  results  to  be  used  in  other  radiative  transfer  applications  such 
as  for  solar  loading  and  energy  balance  calculations. 
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An  updated  version  of  the  AGAUS  Mie  scattering  routine  was  also  used  in  the 
calculations.  This  model  was  introduced  to  improve  the  accuracy  of  forward 
scattered  radiation  predictions  for  large  size  parameter  aerosols. 

Documentary  improvements  include  additional  tables  detailing  the  net  aerosol 
densities  used  in  the  phase  function  calculations,  tables  of  the  real  and  imaginary 
indices  of  refraction  used  in  the  calculations,  and  an  updated  snow  particle  size 
distribution. 

Conclusions 

The  PFNDAT  database  is  an  accurate  description  of  aerosol  scattering  properties 
of  a  wide  range  of  particulates  to  be  encountered  on  the  dirty  battlefield.  This 
upgraded  version  provides  significantly  more  resolution  and  detail  concerning 
these  aerosols  of  interest  to  Army  systems  developers  and  radiative  transfer 
specialists  than  previous  editions. 


1.  INTRODUCTION 

1.1  PFNDAT  Overview 

The  propagation  of  electro-magnetic  energy  within  the  Earth’s  atmosphere 
depends  on  the  wavelength  of  the  radiation  and  on  the  nature  of  the  medium 
being  traversed.  This  medium  consists  of  various  molecular  species  and  aerosol 
particles.  In  this  report  we  describe  the  nature  of  a  particular  set  of  aerosol 
scatterer  classes  related  to  near-ground  propagation  issues  relevant  to  the  Army. 

To  characterize  an  aerosol  species  we  must  be  able  to  identify  the  number,  size, 
shape,  composition,  and  distribution  of  the  aerosol  particles.  In  general  a  given 
aerosol  species  will  consist  of  a  statistical  distribution  of  particle  sizes  and  mean 
real  and  imaginary  indices  of  refraction.  Some  species  are  considered  to  be 
composed  of  weighted  sums  of  more  than  one  particle  type,  as  in  the  case  of 
various  dusts. 

Once  the  nature  of  the  scatterer  particles  is  known,  a  model  can  be  used  to 
determine  the  effect  of  this  species  on  atmospheric  propagation.  However, 
in  general  it  is  not  realistic  to  rely  on  the  direct  properties  of  the  scattering 
species  in  most  radiative  transfer  models.  Instead,  an  intermediate  procedure 
is  used.  This  procedure  determines  the  overall  scattering  properties  of  each 
class  of  scatterer.  Since  each  atmospheric  constituent  scatters  or  absorbs  the 
incident  radiation  according  to  its  own  properties,  and  the  incident  radiation 
may  have  been  previously  scattered  by  another  constituent,  in  order  to  make 
the  atmospheric  propagation  problem  tractable,  a  preprocessing  step  is  often 
needed  where  the  single  scattering  properties  of  a  class  are  determined.  These 
single  scattering  results  can  then  be  used  in  the  more  general  code  to  determine 
multiple  scattering  problems. 

To  evaluate  the  single  scattering  properties  of  a  given  scattering  species  a  number 
of  simplifying  assumptions  axe  often  made: 

•  The  particles  are  assumed  spaced  fax  enough  apart  that  radiation  scattered 
by  one  particle  does  not  affect  how  radiation  is  scattered  from  another 
particle.  Each  scattering  event  is  therefore  independent. 

•  A  Mie  scattering  code  is  assumed  valid  for  predicting  the  behavior  of  each 
scattering  species.  The  usage  of  a  Mie  scattering  method  assumes  that  the 
particles  can  be  approximated  as  spherical  in  shape.  This  assumption  is  often 
made  even  for  nonspherical  particles,  since  in  most  cases  the  orientation 
of  the  particles  is  random  because  no  external  influences  such  as  strong 
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magnetic  fields  or  hydrodynamic  forces  are  present.  Orientation  averaging 
then  produces  nearly  the  same  result  as  the  assumption  of  sphericity. 

•  The  scattering  properties  of  a  given  type  of  particulate  distribution  can  be 
represented  by  a  weighted  integral  over  the  particle  size  distribution.  This 
approach  is  dependent  on  the  assumption  of  independent  scattering  above. 

•  The  particulate  size  distribution  divided  by  the  density  (though  not 
necessarily  the  density  itself)  is  constant  over  the  volume  concerned. 

Having  preprocessed  the  scattering  properties  of  various  aerosol  species,  on  can 
then  accurately  describe  the  scattering  and  absorption  of  radiation  of  a  given 
wavelength  as  it  passes  through  the  atmosphere.  The  relevant  information 
needed  to  determine  these  radiation  results  includes  the  angular  scattering 
probability  distribution,  the  volume  extinction  coefficient  (^),  and  the  single 
scattering  albedo  (w).  The  volume  extinction  coefficient  yd  that  determines  the 
attenuation  of  the  incident  radiation,  is  composed  of  two  parts: 

•  a  scattering  coefficient  /?s  that  describes  the  radiation  scattered  out  of  the 
line  of  sight  (LOS)  without  a  change  in  wavelength 

•  the  absorption  coefficient  yda  that  describes  the  amount  of  radiation  along 
the  LOS  converted  into  other  forms  of  energy  or  that  undergoes  a  change  in 
wavelength. 

These  two  quantities  are  related  to  ^  and  w  by 

W  =  PaUPs  +  /da),  /da  =  IXTyd,  yda  =  (1  -  tu)^.  (1) 

w  represents  the  probability  that  interacting  radiation  will  be  scattered  rather 
than  absorbed:  for  pure  scattering,  =  1;  for  total  absorption,  w  =  0.0. 
The  angular  scattering  distribution  (the  phase  function)  gives  the  directional 
distribution  of  radiation  scattered  by  the  aerosol  under  consideration:  the  phase 
function  P  is  proportional  to  the  probability  that  incoming  radiation  that 
scatters  is  scattered  through  an  angle  6  into  an  element  of  solid  angle  dQ.  The 
phase  function  for  incident  unpolarized  radiation  used  here  is  normalized  as, 

—  /  Pie)  dft  =  -^  /  d4>  f  P{e)  sin(d)  de 

47r  Att  47r  Jo  Jo 

P{e)  sin(d)  d9 

=  1.0,  (2) 


where  8  is  the  scattering  angle. 


1.2  Availability 

EOSAEL  is  available  to  U.S.  Government  Agencies,  specified  allied  organizations, 
and  their  authorized  contractors  at  no  cost.  U.S.  Government  agencies  needing 
EOSAEL  should  send  a  letter  of  request,  signed  by  a  branch  chief  or  division 
director,  to  the  U.S.  Army  Research  Laboratory  (ARL).  Contractors  should  have 
their  Government  contract  monitor  send  the  letter  of  request.  Allied  nation 
organizations  must  request  EOSAEL  through  their  national  representative. 
Please  include,  within  security  restrictions,  a  short  description  of  your  intended 
use(s). 

Release  of  EOSAEL  requires  a  Memorandum  of  Agreement  (MOA)  between 
ARL  and  the  recipient’s  organization.  We  will  send  an  MOA  to  you  for  signature, 
when  you  return  that  to  us  we  sign  it  and  return  a  copy  of  the  MOA  to 
you.  EOSAEL  is  currently  distributed  through  the  DoD  TECNET  facility;  this 
Test  and  Evaluation  Community  NETwork  system  is  located  at  Patuxent  River 
Maryland.  If  you  do  not  already  have  an  account  on  TECNET  we  will  sponsor 
an  account  for  you  and  include  an  application  for  you  to  fill  out.  Return  the 
application  to  ARL  and  we  will  complete  the  account  application  process  for  you. 
You  will  receive  information  about  how  to  log  onto  the  TECNET  (through  the 
Internet,  or  dial-up)  directly  from  TECNET.  If  you  need  additional  help  locating 
or  downloading  EOSAEL  files  after  you  get  your  account,  contact  ARL. 

On  TECNET,  the  EOSAEL  source  code,  DOS  executables,  sample  input  and 
output  files  are  available.  Documentation  for  the  modules  is  included  as 
postscript  files  suitable  for  viewing  or  printing. 

Specific  technical  questions  concerning  PFNDATshould  be  directed  to  David 
Tofsted  at  U.S.  Army  Research  Laboratory,  (505)  678-3039  commercial  and 
258-3039  DSN,  or  via  e-mail  at  dtofsted@arl.mil. 

The  EOSAEL  point  of  contact  at  ARL  is  Dr.  Alan  Wetmore. 

1.2.1  Mailing  Address 

U.S.  Army  Research  Laboratory 
ATTN:  AMSRL-IS-ES  (Dr.  Wetmore) 

2800  Powder  Mill  Road 
Adelphi,  MD  20783-1197 

1.2.2  Phone  and  Electronic  Mail 

COM:(301)  394-2499 
FAX:(301)  394-4797 
DSN:  290-2499 

emaihawetmoreOarl  .mil 
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2.  TECHNICAL  DOCUMENTATION 

2.1  EOSAEL  Application  -  PFNDAT 

The  Electro- Optical  Systems  Atmospheric  Effects  Library  (EOSAEL)  is  a 
coordinated  set  of  databases  and  models  that  allow  for  characterization  of 
numerous  atmospheric  effects  on  electro-optical  systems.  Because  of  the 
pervasive  nature  of  aerosol  scattering  effects  in  all  aspects  of  propagation 
problems,  and  due  to  the  desire  to  not  have  to  compute  atmospheric  scattering 
properties  on  a  case-by-case  basis,  it  is  reasonable  to  accumulate  a  collection 
of  scattering  phase  functions  for  a  wide  variety  of  atmospheric  conditions  in 
one  place.  The  collection  is  called  PFNDAT,  the  Phase  FunctioN  DATabase. 
PFNDAT  contains  many  commonly  encountered  aerosol  types  within  the 
atmospheric  boundary  layer,  including  the  main  haze  aerosols  (rural,  urban, 
and  maritime),  two  classes  of  fog  (radiative  and  advective),  as  well  as 
precipitation  classes  (drizzle,  rain,  and  one  type  of  snow)  and  battlefield  induced 
contaminant  aerosols  (fog  oil,  hexachloroethane  (HC),  white  phosphorous  (WP), 
and  dusts).  This  set  covers  the  scattering  effects  of  the  majority  of  cases  of 
aerosols  encountered  by  the  Army  within  the  near  surface  environment,  up  to 
approximately  1000  m  above  ground  level.  Various  scattering  models  for  upper 
air  aerosols  (most  cloud  types)  are  not  included  in  this  database;  databases 
developed  by  the  Air  Force  address  this  area. 

As  a  result  of  the  development  of  this  common  database,  many  models  utilize 
PFNDAT  as  model  input.  Some  EOSAEL  modules  use  the  volume  extinction 
coefficients  to  determine  transmission  along  lines  of  sight;  two  EOSAEL  modules 
(the  Approximate  Multiple  Scattering  Module  (ASCAT)  and  the  Finite  CLOUD 
contrast  transmission  module  (FCLOUD))  deal  with  aspects  of  the  scattered 
radiation  and  therefore  use  the  phase  function  information.  The  Weather  and 
Atmospheric  Visualization  Effects  for  Simulation  suite  of  codes  accesses  the 
PFNDAT  phase  function  and  extinction  information  when  computing  scattering 
effects  in  the  Boundary  Layer  Illumination  and  Radiation  Balance  radiative 
transfer  model. 

Scattering  results  for  38  different  aerosol  distributions  have  been  included  in 
the  data  base  to  cover  the  aerosol  environments  expected  for  Army  purposes. 
The  extinction  coefficients  and  phase  functions  produced  were  generated  with 
the  computer  code  AGAUS  (Miller  1983),  that  uses  the  classical  Mie  scattering 
theory  approach. 

Phase  functions  and  extinction  coefficients  for  the  fogs,  rains,  snow,  and 
maritime,  urban  and  rural  aerosols  were  generated  at  wavelengths  of  0.35,  0.40, 
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0.45,  0.50,  0.55,  0.60,  0.65,  0.70,  0.75,  1.06,  3.0,  3.5,  4.0,  4.5,  5.0,  8.0,  8.5,  9.0, 

9.5,  10.0,  10.5,  11.0,  11.5,  12.0,  14.0,  16.0,  18.0,  20.0,  25.0,  30.0,  35.0,  and  40.0 
/j,m.  The  manmade  aerosols  were  generated  at  wavelengths  of  0.55,  1.06,  3.0, 

3.5,  4.0,  4.5,  5.0,  8.0,  8.5,  9.0,  9.5,  10.0,  10.5, 11.0,  11.5,  and  12.0  fim.  The  phase 
functions  were  computed  at  65  angles  spaced  unequally  between  0.0  and  180.0 
degrees.  The  selection  of  output  angles  for  each  aerosol  case  provides  more  values 
of  the  phase  function  in  which  the  variation  is  most  rapid,  particularly  in  the 
forward  and  backward  peak  directions  (0°  and  180°,  respectively).  The  phase 
functions  for  all  of  the  distributions  described  here  are  presented  graphically  in 
appendix  B. 

2.1.1  Phase  Function  Production 

The  following  symbology  and  terms  are  used  in  the  remainder  of  this  report.  The 
calculation  of  phase  function  information  for  natural  aerosols  usually  entails  a 
combination  of  effects  from  more  than  one  scattering  substance  type,  for  which 
each  type  may  have  its  own  relative  density  distribution  and  refractive  index 
properties.  For  example,  dusts  are  composed  of  various  component  quantities 
of  quartz,  montmorillonite,  and  ammonium  sulfate  particles.  Each  particle  type 
has  a  particulate  bulk  density  p  (in  units  of  g/cm®)  for  a  given  volume  of  the 
scattering  substance.  However,  each  species  must  also  be  characterized  by  a 
given  mass  concentration  C,  also  given  in  units  of  density  (g/cm®),  representing 
the  weight  of  lofted  material  mixed  within  a  unit  volume  of  air.  C  is  often 
referred  to  by  the  term  liquid  water  content  from  a  meteorological  standpoint. 
The  term  mass  concentration  is  used  to  describe  both  effects  in  this  text. 

Relating  the  quantities  p  and  C,  the  particle  size  distribution  is  denoted  by  n(r), 
where  r  is  the  radius  of  a  particle  (using  Mie  theory  all  particles  are  considered 
spherical).  The  total  number  of  particles  per  unit  volume  is  denoted  by  the 
number  density  N,  which  relates  to  the  particle  size  distribution  through 

r°°  dN 

N=  n(r)  dr;  n(r)  =  .  (3) 

Jo 

n(r)  has  units  of  particles  per  cm^-pm. 

Writing  the  mass  concentration  as  a  function  of  the  particle  size  distribution, 

C  =  y  n(r)  p  dr,  (4) 

where  (A/Z)'Kr^  p  is  the  mass  of  a  particle  of  radius  r. 

The  particle  size  distributions  are  comprised  of  the  Air  Force  Geophysics 
Laboratory  (AFGL)  maritime,  urban,  and  rural  models  at  relative  humidities 
of  0,  50,  70,  80,  90,  95,  98,  and  99  percent,  the  AFGL  fog  models  for  heavy 
advection  and  radiation  types  (Shettle  and_  Fenn  1979);  three  rain  models 
(drizzle,  widespread,  and  thunderstorm),  a  snow  model,  moderate  and  heavy 
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aerosol  dust  models,  a  high-explosive  (HE)  dust  model,  and  three  smoke  types. 
The  WP  smoke  type  is  calculated  at  three  relative  humidities  (17,  50,  and  90 
percent).  The  rain  models  use  the  Marshall-Palmer  (MP)  distribution  based 
on  the  work  of  Marshall  and  Palmer  (1948).  The  snow  and  fog  models  use  the 
modified  gamma  (MG)  distribution.  The  smoke  distributions  are  lognormal. 
The  dust,  rural,  urban,  and  maritime  distributions  are  bimodal  lognormal  (large 
and  small  particle  component  lognormal  distributions). 

2.1.2  Particle  Size  Distribution  Models  Used 


The  following  aerosol  particle  size  distributions  are  used  in  evaluating  the  phase 
functions  contained  in  the  PFNDAT  database. 

Marshall-Palmer  Rain  Distribution 

The  MP  distribution  is  given  in  the  AGAUS  documentation  (The  1982  AGAUS 
documentation  lists  B  incorrectly  as  8.2  x  10“^,  while  the  value  in  equation  (5) 
is  consistent  with  Pruppacher  and  Klett  (1980).)  by  the  equation 


/  N  <1-^  ..  /  7-,  N 

”(*')  =  ^  exp(— Hr) 

dr 

(5) 

A  =  1.6  X  10"®  cm"®^m"^, 

(6) 

H  =  8.2  X  10-®R-°-2Vm"\ 

(7) 

where  R  is  the  rain  rate  in  mm/hr.  Using  this  form  and  equation  (3)  yields  the 
relationship, 

N  =  A/B  =  1951  X  10"®  cm’^  (8) 

The  mass  concentration  for  this  distribution  is  given  by, 

=  8  TT  p  ^  =  89  X  10-=*  i^®■«^g/m^  (9) 


where  p  =  1  g/cm®  is  used  for  the  bulk  density  of  water. 

Modified  Gamma  Distribution 


The  modified  gamma  (MG)  distribution  is  given  in  Shettle  and  Fenn  (1979)  by 
the  equation 


dN 

—  =  A  r°‘  exp  [— 6r*’'] , 
dr 


(10) 


where  A,  b,  a  and  7  are  fit  coefficients  of  the  distribution.  A  has  units  of 
number  of  particles  per  cubic  centimeter  per  /zm""*"*.  A  similar  equation  is  used 
by  AGAUS: 


n(r) 


dN 

dr 


=  Vcr 


(11) 


where  b  must  equal  a/{‘yr2)  to  convert  from  the  form  of  the  exponential 
argument  in  equation  (10)  to  the  form  in  equation  (11). 
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Equation  (11)  is  the  distribution  contained  within  the  AGAUS  program.  The 
equation  in  the  1983  AGAUS  documentation  incorrectly  contained  an  r  instead 
of  a  7  in  the  denominator  of  the  first  term  of  the  exponential  argument.  However 
this  equation  is  in  error,  since  it  does  not  contain  any  particle  density  dependence 
as  in  Pruppacher  and  Klett  (1980).  A  direct  integration  of  the  particle  size 
distribution,  using  definition  3.478.1  from  Gradshteyn  and  Ryzhik  (1980),  allows 
us  to  substitute  a  function  for  A,  which  depends  on  6,  a,  and  7.  The  resulting 
equation  matches  that  listed  in  Pruppacher  and  Klett  (1980): 


(12) 

(13) 

(14) 


Solving  this  equation  for  A  and  substituting  the  function  in  Tc,  a,  and  7  for  5, 
a  new  equation  for  the  particle  size  distribution  equation  is  obtained. 


where. 


(15) 


(16) 


The  mass  concentration  for  this  distribution  can  be  expressed  by 


C  = 


ri»>  +  (3/7)1 

/  \  3/7 

(?)  r(tf) 


(17) 


Lognormal  Distribution 

The  lognormal  distribution  is  given  in  Shettle  and  Fenn  (1979)  by  the  equation 


dN  _ 

N 

exp 

(logior-logiorg)2' 

dr 

.ln(10)r  V^asF . 

2  (t|p 

(18) 


N  is  the  aerosol  particle  number  density  (particles  per  cm^),  Vg  is  the  distribution 
geometric  mean  radius  (or  mode  radius)  in  fiia,  and  usp  Is  the  width  of  the 
distribution  measured  in  logio  space. 

A  similar  equation  is  used  for  the  lognormal  distribution  in  the  AGAUS  program, 
but  with  a  different  meaning  for  the  a  term.  In  AGAUS 


dN 

N 

exp 

(In  r  —  In  Tj  'f‘ 

dr 

r  In(crj) 

(19) 


where  ag  is  called  the  geometric  mean  standard  deviation.  The  relationship 
between  ag  and  asp  is 

(TsF  =  logio(<73);  CTSF  In(lO)  =  In(aj).  (20) 

The  above  lognormal  distribution  equation  is  correct  for  the  code  contained 
within  the  AGAUS  program.  (The  equation  contained  in  the  1983  AGAUS 
documentation  was  missing  the  appropriate  normalization  terms.  Perhaps  this 
is  intentional,  since  AGAUS  appears  to  use  a  different  method  of  normalization.) 
The  mass  concentration  equation  for  this  distribution  is  given  by 


2.2  Maritime,  Urban,  and  Rural  Aerosol  Models 

The  maritime,  urban,  and  rural  aerosol  models  are  identical  to  those  found 
in  Shettle  and  Fenn  (1979)  and  are  bimodal  lognormal,  with  the  mode  radius 
varying  as  a  function  of  relative  humidity.  The  rural  aerosol  model  consists 
of  small  and  large  rural  distributions  with  correspondingly  different  indices  of 
refraction.  Similarily,  the  urban  aerosol  model  consists  of  small  and  large  urban 
distributions.  The  maritime  aerosol  model  consists  of  the  small  rural  distribution 
along  with  a  large  particle  continental  oceanic  distribution. 

The  indices  of  refraction  for  the  individual  aerosols  are  in  Appendix  A. 

The  number  densities  (rounded)  for  each  mode  of  the  distribution  type  along 
with  the  mode  radius  and  variance  data,  the  resulting  extinction  coefficient  data, 
and  the  related  liquid  water  content  information  are  provided  for  the  user  who 
may  wish  to  change  the  visibilities.  The  number  density,  mode  radius  (vg),  and 
variance  {ag)  information  for  the  maritime  cierosols  are  contained  in  table  1 
and  the  extinction  data  is  contained  in  table  2.  The  same  information  for  the 
urban  aerosols  is  in  tables  3  and  4.  The  same  information  for  the  rural  aerosols  is 
contained  in  tables  5  and  6.  The  liquid- water-content  (mass  concentration)  data 
is  contained  in  table  7.  The  latter  information  is  derived  using  equation  (21) 
for  each  aerosol  type.  The  extinction  coefficients  in  these  tables  were  generated 
for  a  5.0-km  visibility,  with  number  densities  corresponding  to  those  found  in 
tables  8  through  10  of  Shettle  and  Fenn  (1979). 
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Table  1.  Mode  radii  (/um),  spread,  and  number  densities  (cm“^)  as  functions  of 
relative  humidity  for  the  small  (S)  and  large  (L)  modes  of  the  AFGL  maritime 
haze  aerosol  model. 


Relative  Humiditvf%)  1 

Qty 

0 

50 

70 

80 

90 

95 

98 

99 

NiS) 

38251 

35129 

27757 

13902 

9697 

6976 

4360 

2948 

N(L) 

386.4 

354.8 

280.4 

140.4 

98.0 

70.5 

44.0 

29.8 

rg  (S) 

0.02700 

0.02748 

0.02846 

0.03274 

0.03884 

0.04238 

0.04751 

0.05215 

r,  (L) 

0.1600 

0.1711 

0.2041 

0.3180 

0.3803 

0.4606 

0.6024 

0.7505 

(S) 

2.239 

2.239 

2.239 

2.239 

2.239 

2.239 

2.239 

2.239 

(L) 

2.512 

2.512 

2.512 

2.512 

2.512 

2.512 

2.512 

2.512 

Table  2.  Extinction  coefficients  (km  versus  wavelength  (//m)  for  the  maritime 
aerosol  model 


A 

Relative  Humiditvf%)  1 

(fim) 

0 

50 

70 

80 

90 

95 

98 

99 

0.35 

0.9605 

0.9499 

0.9163 

0.8549 

0.8582 

0.8412 

0.8217 

0.8117 

0.40 

0,8991 

0.8953 

0.8700 

0.8257 

0.8325 

0.8200 

0.8045 

0.7970 

0.45 

0.8530 

0.8478 

0.8341 

0.8057 

0.8078 

0.8018 

0.7925 

0.7872 

0.50 

0.8083 

0.8076 

0.7969 

0.7849 

0.7874 

0.7861 

0.7797 

0.7780 

0.55 

0.7696 

0.7718 

0.7711 

0.7691 

0,7714 

0,7721 

0.7695 

0.7694 

0.60 

0.7324 

0.7383 

0.7439 

0.7581 

0.7566 

0.7569 

0.7565 

0.7600 

0.65 

0.7049 

0.7105 

0.7201 

0.7414 

0.7425 

0.7494 

0.7540 

0.7576 

0.70 

0,6768 

0.6867 

0.6977 

0.7285 

0,7306 

0.7397 

0.7531 

0.7571 

0.75 

0.6559 

0.6596 

0.6743 

0.7208 

0.7216 

0.7301 

0.7398 

0.7460 

1.06 

0.5430 

0.5566 

0.5946 

0,6738 

0.6783 

0.6995 

0.7198 

0.7321 

3.00 

0.3265 

0.3414 

0.3900 

0.5236 

0.5541 

0.5992 

0.6563 

0.6881 

3.50 

0.2410 

0.2618 

0.3234 

0.4937 

0.5325 

0.5975 

0.6742 

0.7248 

4.00 

0.2119 

0.2273 

0.2747 

0.4271 

0.4700 

0.5397 

0.6301 

0.6943 

4.50 

0.1917 

0.2040 

0.2440 

0.3815 

0.4242 

0.4950 

0.5902 

0.6621 

5.00 

0.1631 

0.1747 

0.2119 

0.3436 

0.3879 

0.4594 

0.5598 

0,6368 

8.00 

0.07913 

0.08616 

0.1096 

0.2003 

0.2377 

0.2974 

0.3930 

0.4768 

8.50 

0.09768 

0.1018 

0.1177 

0.1902 

0.2215 

0.2751 

0.3646 

0.4464 

9.00 

0.1207 

0.1229 

0.1317 

0.1844 

0.2099 

0.2564 

0.3369 

0.4147 

9.50 

0.09504 

0.09732 

i  0.1067 

0.1586 

0.1828 

0.2266 

0.3023 

0.3761 

10.00 

0.08018 

0.08099 

0.08988 

0.1371 

0.1592 

0.1981 

0.2671 

0.3350 

10.50 

0.06671 

0.06866 

0.07741 

0.1206 

0.1403 

!  0.1741 

0.2340 

0,2935 

11.00 

0.05683 

0.05982 

0.07113 

0,1186 

0.1379 

0,1704 

0.2240 

0.2767 

11.50 

0.05261 

0.05690 

0.07234 

0.1295 

0.1981 

0,1848 

0.2382 

0.2878 

12.00 

0.04529 

0.05268 

0.07668 

0.1519 

0.2072 

0.2165 

0.2747 

0.3276 

14.00 

0.03274 

0,04895 

0.09399 

0.2116 

0.2488 

0.2992 

0.3717 

0.4335 

15.00 

0.03824 

0.05470 

0.1002 

0.2208 

0.2596 

0.3126 

0.3881 

0.4526 

18.00 

0.05800 

0.07020 

0.1070 

0.2184 

0.2571 

0.3125 

0.3944 

0.4654 

20.00 

0.05147 

0.06235 

0.09504 

0.1968 

0.2338 

0.2882 

0.3709 

0.4448 

25.00 

0.04225 

0.05017 

0.07445 

0.1557 

0.1880 

0.2369 

0.3152 

0.3886 

30.00 

0.04105 

0.04662 

0.06401 

0.1265 

0.1534 

0.1956 

0.2667 

0.3361 

35.00 

0.04729 

0.05128 

0.06379 

0.1120 

0.1336 

0.1696 

0.2321 

0.2950 

40.00 

0.06735 

0.07085 

0.07976 

0,1150 

0.1320 

0.1631 

0.2195 

0.2773 

Number 

density 
(cm  ) 

38637 

35484 

28037 

14042 

9795 

7047 

4404 

2978 

Table  3.  Mode  radii  spread,  and  number  densities  (cm“^)  as  functions 

of  relative  humidity  for  the  small  (S)  and  large  (L)  modes  of  the  AFGL  urban 
haze  aerosol  model. 


Relative  Humiditvf%l  1 

Qty 

0 

50 

70 

80 

90 

99 

Ar(S) 

87204 

83354 

64829 

42776 

27693 

18217 

7286 

A'(L) 

10.9 

10.4 

8.1 

5.4 

3.5 

2.3 

1.3 

(S) 

0.02500 

0.02563 

0.02911 

0.03514 

0.05996 

(L) 

0.4000 

0.4113 

0.4777 

0.5805 

1.1690 

cTff  (S) 

2.239 

2.239 

2.239 

2.239 

2.239 

2.239 

2.239 

2.239 

(L) 

2.512 

2.512 

2.512 

2.512 

2.512 

2.512 

2.512 

2.512 

Table  4.  Extinction  coefficients  (km  versus  wavelength  (/um)  for  the  urban 
aerosol  model 


A 

Relative  Humiditv(%)  1 

(^m) 

0 

50 

70 

80 

90 

95 

98 

99 

0.35 

1.141 

1.144 

1.149 

1.144 

1.121 

1.083 

1.024 

0.9800 

0.40 

1.031 

1.032 

1.036 

1.034 

1.018 

0.9941 

0.9547 

0.9254 

0.45 

0.9334 

0.9345 

0.9358 

0.9353 

0.9267 

0.9133 

0.8900 

0.8724 

0.50 

0.8473 

0.8481 

0.8478 

0.8489 

0.8458 

0.8403 

0.8280 

0.8212 

0.55 

0.7711 

0,7717 

0.7714 

0.7723 

0,7714 

0.7717 

0.7694 

0.7704 

0.60 

0.7040 

0.7045 

0.7029 

0.7046 

0.7082 

0.7112 

0.7181 

0.7233 

0.65 

0.6453 

0.6457 

0.6441 

0,6463 

0.6519 

0.6581 

0.6707 

0.6810 

0.70 

0,5944 

0,5946 

0.5929 

0.5951 

0.6004 

0.6089 

0.6259 

0,6398 

0.75 

0.5473 

0.5474 

0.5459 

0.5474 

0.5545 

0.5645 

0.5847 

0.6016 

1.06 

0.3632 

0.3551 

0.3596 

0.3571 

0.3618 

0.3737 

0.3982 

0.4237 

3.00 

0.1244 

0.1311 

0.1651 

0.2046 

0.2317 

0.2508 

0.2727 

0.2909 

3.50 

0.1126 

0.1130 

0.1147 

0.1136 

0.1133 

0.1168 

0.1278 

0.1433 

4.00 

0.1042 

0.1039 

0.1040 

0.10053 

0.09799 

0.09921 

0.1069 

0.1196 

4.50 

0.09813 

0.09782 

0.09804 

0.09489 

0.09248 

0.09359 

0.1009 

0.1130 

5.00 

0.09165 

0.09139 

0.09138 

0.08823 

0,08601 

0.08713 

0.09420 

0.1061 

8.00 

0.06710 

0,06670 

0,06597 

0.06451 

0.06443 

0.06737 

0.07621 

0.08904 

8.50 

0.08252 

0.08175 

0.07797 

0.07135 

0.06749 

0.06781 

0.07483 

0.08690 

9.00 

0.09358 

0.09427 

0.09569 

0.08812 

0.07888 

0.07437 

0.07643 

0.08641 

9.50 

0.08207 

0.08212 

0.08108 

0.07400 

0.06752 

0.06531 

0.06965 

0.08059 

10.00 

0.07552 

0.07544 

0.07401 

0.06764 

0.06196 

0.05818 

0.06494 

0.07559 

10.50 

0.07069 

0.07056 

0,06899 

0.06305 

0.05807 

0.05637 

0.06722 

0.07097 

11.00 

0.06619 

0.06616 

0.06544 

0.06122 

0.05775 

0.05719 

0.06193 

0.07144 

11.50 

0.06310 

0.06328 

0.06388 

0.06224 

0.06110 

0.06215 

0.06790 

0.07733 

12.00 

0.06061 

0.06114 

0.06405 

0.06644 

0.06843 

0,07007 

0.07918 

0,08958 

14.00 

0.05348 

0.05523 

0.06538 

0.07757 

0.08652 

0.09431 

0.1058 

0.1186 

15.00 

0.05554 

0.05754 

0.06845 

0.08050 

0.08904 

0.09641 

0.1078 

0.1207 

18.00 

0.05115 

0.05285 

0.06236 

0.07300 

0.08085 

0.08786 

0.09962 

0.1132 

20.00 

0.04961 

0.05102 

0.05903 

0.06741 

0.07356 

0.07978 

0.09083 

0.1042 

25.00 

0.04139 

0.04256 

0.04903 

0.05549 

0.06048 

0.06598 

0.07666 

0.09000 

30.00 

0.03646 

0.03737 

0.04237 

0.04726 

0.05120 

0.05611 

0.06631 

0.07954 

35.00 

0,03348 

0.03431 

0.03874 

0.04279 

0.04604 

0.05034 

0.05984 

0.07255 

40.00 

0.03119 

0.03206 

0,03658 

0.04075 

0.04403 

0.04822 

0.05736 

0.06971 

Number 

density 

87215 

83364 

64837 

42781 

27697 

18219 

10517 

7287 

(cn^ 

• 
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Table  5.  Mode  radii  spread,  and  number  densities  (cm”^)  as  functions  of 

relative  humidity  for  the  small  (S)  and  large  (L)  modes  of  the  AFGL  rural  haze 
aerosol  model. 


Relative  Humiditvf%)  1 

Qty 

0 

50 

70 

80 

90 

95 

98 

99 

N(S) 

79076 

76305 

70804 

51674 

33895 

27052 

19290 

14761 

N{L) 

9.9 

9.5 

8.9 

6.4 

4.2 

3.4 

2.4 

1.9 

rg  (S) 

0.02700 

0.02748 

0.02846 

0.03274 

0.03884 

0.04238 

0.04751 

0.05215 

rg  (L) 

0.4300 

0.4377 

0.4571 

0.5477 

0.6462 

0.7078 

0.9728 

1.1760 

0-3  (S) 

2.239 

2.239 

2.239 

2.239 

2.239 

2.239 

2.239 

2.239 

(^g  (L) 

2.512 

2.512 

2.512 

2.512 

2.512 

2.512 

2.512 

2.512 

Table  6.  Extinction  coefficients  (km  versus  wavelength  (/um)  for  the  rural 
aerosol  model 


A 

Relative  Humiditvf%)  1 

(/xm) 

0 

50 

70 

80 

90 

95 

98 

99 

0.35 

1.201 

1.200 

1.199 

1.187 

1.153 

1.135 

1.091 

1.063 

0.40 

1.071 

1.071 

1.069 

1.062 

1.039 

1.027 

0.9990 

0.9819 

0.45 

0.9571 

0.9565 

0.9548 

0.9534 

0.9398 

0.9314 

0.9153 

0.9071 

0.50 

0.8574 

0.8572 

0.8565 

0.8569 

0.8509 

0.8468 

0,8401 

0.8376 

0.55 

0.7711 

0.7714 

0.7711 

0.7730 

0.7700 

0.7714 

0.7703 

0.7730 

0.60 

0.6958 

0.6959 

0.6962 

0.7009 

0.7021 

0.7040 

0.7092 

0.7158 

0.65 

0.6305 

0.6306 

0.6317 

0.6380 

0.6419 

0.6457 

0.6553 

0,6655 

0.70 

0.5740 

0.5743 

0.5755 

0.5830 

0.5871 

0,5922 

0.6077 

0.6204 

0.75 

0.5216 

0.5220 

0.5235 

0.5327 

0.5374 

0.5450 

0.5631 

0.5781 

1.06 

0.3236 

0.3235 

0.3259 

0.3341 

0.3400 

0.3473 

0.3713 

0.3936 

3.00 

0.08639 

0.09261 

0.1069 

0.1565 

0.1957 

0.2127 

0.2711 

0.2738 

3.50 

0.08051 

0.08100 

0.08447 

0.09481 

0.09666 

0.09921 

0,1290 

0.1480 

4.00 

0.07271 

0.07427 

0.07694 

0.08454 

0.08379 

0.08501 

0.1140 

0.1321 

4.50 

0.07091 

0.07096 

0.07344 

0.08045 

0,07970 

0.08070 

0.1096 

0.1277 

5.00 

0.06524 

0.06537 

0.06929 

0.07443 

0.07384 

0,1060 

0.1278 

0.1418 

8.00 

0.03321 

0.03382 

0.03615 

0.04611 

0.05052 

0.05324 

0.08235 

0,1018 

8.50 

0.06738 

0.06641 

0.06620 

0.06550 

0.06093 

0.06045 

0.08427 

0.1015 

9.00 

0.09769 

0.09736 

0.09862 

0.09676 

0.08469 

0.07981 

0.09496 

0.1074 

9.50 

0.08006 

0.07931 

0.07954 

0.07641 

0.06736 

0,06457 

0.08152 

0.09548 

10.00 

0.07066 

0.06982 

0.06999 

0.06711 

0.05972 

0.05765 

0.07397 

0,08749 

10.50 

0.06362 

0.06290 

0.06293 

1  0.06039 

0.05439 

0.05283 

0.06768 

0,08015 

11.00 

0.05755 

0.05697 

0.05715 

0.05583 

0.05204 

0.05148 

0.06587 

0.07790 

11.50 

0.05311 

0.05275 

0.05330 

0.05415 

0.05326 

0.05414 

0.06959 

0.08203 

12.00 

0.04982 

0.04979 

0.05107 

0.05560 

0.05860 

0.06090 

0.07946 

0.09342 

14.00 

0.04109 

0.04244 

0.04645 

0.06193 

0.07334 

0.07887 

0.1049 

0.1226 

15.00 

0.04517 

0.04684 

0.05142 

0.06780 

0.07804 

0.08297 

0.1091 

0,1269 

18.00 

0.04798 

0.04882 

0.05214 

0.06436 

0.07189 

0.07591 

0.1028 

0,1214 

20.00 

0.04854 

0.04918 

0.05191 

0.06168 

0.06675 

0.06980 

0.09562 

0.1138 

25.00 

0.03996 

0.04057 

0.04289 

0.05116 

0.05500 

0.05744 

0.08165 

0.09920 

30.00 

0.03542 

0.03586 

0.03769 

0,04416 

0.04688 

0.04882 

0.07092 

0.08751 

35.00 

0.03422 

0.03452 

0,03602 

0.04115 

0.04292 

0.04436 

0.06406 

0.07938 

40.00 

0.03320 

0.03352 

0.03499 

0.03985 

0.04150 

0.04286 

0.06140 

0.07583 

Number 

density 

79086 

76315 

70813 

51680 

33899 

27055 

19292 

14763 

Table  7.  Liquid  water  content  as  functions  of  relative  humidity  for  small 

mode,  large  mode,  and  total  content  for  rural,  urban,  and  maritime  aerosols. 


Aerosol  Type 

Relative  Humiditvf%l  1 

0 

50 

70 

80 

90 

95 

99  1 

Maritime 

1 

Small 

58.7 

56.8 

49.9 

38.0 

44.3 

41.4 

36.4 

32.6 

Large 

301.7 

338.7 

454.4 

860.5 

1027.3 

1313.0 

1833.2 

2400.9 

Total 

360.3 

395.5 

504.3 

898.5 

1071.6 

1354.4 

1869.7 

2433.5 

Urban 

1 

Small 

109.4 

124.6 

144.7 

158.4 

167.5 

176.7 

182.3 

Large 

137.9 

168.3 

234.8 

282.1 

395.8 

561.7 

Total 

247.3 

292.9 

393.3 

449.6 

572.5 

744.0 

Rural 

1 

Small 

121.3 

123.4 

127.2 

141.3 

154.8 

160.5 

161.2 

163.2 

Large 

150.0 

151.8 

162.0 

216.0 

229.8 

421.1 

588.9 

Total 

271.3 

275.3 

289.2 

341.7 

370.8 

390.3 

582.3 

752.1 

2.3  Fog  Models 

The  Shettle  and  Fenn  (1979)  data  were  also  used  for  the  heavy  advection 
and  radiation  fog  models.  These  models  use  the  MG  distribution,  and  their 
extinction  coefficients  are  given  in  table  8  for  wavelengths  of  0.35  to  40.0  fim. 
For  heavy  (advection)  fog,  the  mode  radius  was  10.0  /xm,  with  a  number  density 
of  20  particles/cm^,  a  was  set  to  3,  and  7  was  set  to  1;  for  radiation  fog,  the 
mode  radius  was  2.0  /xm,  with  a  number  density  of  200  parti cles/cm*,  7  was 
again  set  to  1,  and  a  was  set  to  6.  The  equivalent  mass  concentrations  for  these 
two  distributions  are  in  table  9  along  with  the  rain  and  snow  results  from  the 
following  section.  Other  details  may  be  found  in  Shettle  and  Fenn  (1979). 

2.4  Rain  and  Snow  Models 

Particle  sizes  of  rain  and  snow  generally  are  quite  large  compared  to  visible 
and  infrared  wavelengths,  making  Mie  calculations  to  determine  phase  functions 
impractical.  According  to  Hodkinson  and  Greenleaves  (1963),  when  the 
airborne  particles  of  an  aerosol  species  are  larger  than  a  few  wavelengths  of 
the  radiation  being  transmitted  and  a  range  of  particle  sizes  or  wavelengths 
exists,  the  combined  single-scattering  characteristics  may  be  approximated 
by  a  combination  of  Fraunhofer  diffraction  and  geometrical  transmission  and 
reflection.  While  studying  light  scattering  by  irregular  particles  larger  than 
the  wavelength  (such  as  snow),  Hodkinson  (1963)  found  that,  although  the 
diffraction  patterns  of  individual  irregular  particles  vary  greatly  with  shape, 
the  resultant  forward  diffraction  lobe  for  an  ensemble  of  nonspherical  particles 
with  random  orientations  would  be  similar  to  an  ensemble  of  spheres  with  cross- 
sectional  areas  equal  to  the  particles’  areas. 
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Table  8.  Extinction  coefficients 


(km  versus  wavelength  (^m)  for  various  fogs 


Wavelength 

Advection 

Radiation 

(^m) 

Fog 

Fog 

0.35 

28.52 

8.404 

0.40 

28.60 

8.478 

0.45 

28.63 

8.533 

0.50 

28.69 

8.580 

0.55 

28.75 

8.614 

0.60 

28.81 

8.690 

0.65 

28.84 

8.714 

0.70 

28.89 

8.799 

0.75 

28.94 

8.840 

1.06 

29.20 

9.108 

3.00 

30.14 

9.608 

3.50 

30.79 

12.739 

4.00 

31.12 

11.925 

4.50 

31.40 

10.292 

5.00 

31.70 

9.049 

8.00 

33.78 

3.992 

8.50 

34.30 

3.400 

9.00 

34.84 

2.884 

9.50 

35.27 

2.453 

10.00 

35.25 

2.093 

10.50 

33.57 

1.919 

11.00 

30.36 

2.081 

11.50 

28.24 

2.498 

12.00 

28.09 

3.073 

14.00 

31.17 

4.743 

15.00 

32.02 

4.965 

18.00 

33.82 

4.822 

20.00 

34.63 

4.199 

25.00 

35.80 

2.989 

30.00 

36.44 

2.160 

35.00 

35.72 

1.765 

40.00 

34.62 

1.673 

Number 

density 

20.00 

200.0 

_ _ 

Table  9.  Liquid  water  content  (/ig/m^)  as  a  function  of  the  type  of  precipitation. 


PrecipitationType 

Liquid  Water  Content  (//g/m^) 

Advection  Fog 

372300 

Radiation  Fog 

15640 

Rain  (Drizzle),!  mm/hr 

89000 

Rain  (Widespread), 5  mm/hr 

344000 

Rain  (Thunderstorms), 10  mm/hr 

615700 

Snow 

68150 

A  special  version  of  AGAUS  for  treating  SNOW  cases  (AGSNOW)  (Deepak  et 
al.  1982)  was  developed  to  calculate  the  single-scattering  characteristics  of 
large  spherical  or  irregular  particles.  AGSNOW  consisted  of  a  combination  of 
Fraunhofer  diffraction,  geometrical  reflection,  and  a  parameterization  (Pollack 
and  Cuzzi  1980)  for  the  refracted  and  internally  reflected  energy.  Mie  theory 
is  used  to  compute  the  phase  functions  for  particles  with  size  parameters  less 
than  a  user  defined  upper  bound,  and  to  compute  extinction  and  scattering 
cross  sections.  The  parameterization  used  for  the  refracted  and  internally 
reflected  energy  precludes  the  calculation  of  effects  of  rainbows  and  other 
such  optical  phenomena;  thus,  the  phase  functions  at  the  affected  angles  are 
only  approximate.  Because  the  phase  functions  for  typical  rain  and  snow 
size  distributions  have  sharp  peaks  in  the  forward  direction,  a  set  of  angles 
concentrated  in  the  forward  direction  must  be  used  for  these  phase  functions. 
Otherwise,  difficulties  arise  in  the  phase  function  interpolation  processes  in 
various  EOSAEL  modules:  a  new  set  of  angles  is  automatically  read  when  the 
rain  and  snow  phase  functions  are  used. 

Phase  functions  for  rain  were  generated  using  the  AGSNOW  code  for  a  MP 
particle  size  distribution  with  rain  rates  of  1,  5,  and  10  mm/h,  corresponding  to 
rain  types  of  drizzle,  widespread,  and  thunderstorm,  respectively,  and  to  number 
densities  of  about  0.0019,  0.0027,  and  0.0031  particles/cm^,  respectively.  The 
MP  distribution  yields  somewhat  larger  particles  at  higher  rain  rates,  causing 
the  forward  direction  lobe  of  the  phase  function  to  be  narrower  and  sharper.  The 
phase  functions  given  are  considered  to  be  reasonable  for  most  rain  conditions. 
The  extinction  coefficient  for  other  rain  rates  may  be  calculated  using  the 
algorithm  given  in  the  EOSAEL  XSCALE  module  (Fiegel  1994). 

Snow  phase  functions  were  calculated  using  an  MG  particle  size  distribution 
that  had  been  fit  to  a  measured  size  distribution  (Unpublished  results  taken 
from  SNOW-ONE-A)  as  shown  in  figure  1.  The  MG  model  was  chosen  because 
it  provides  a  realistic  simulation  of  the  relatively  slow  particle  density  decrease 
for  r  <  Tc.  The  particle  number  density  used  was  N  =  40  particles/m®,  and 
the  peak  density  particle  diameter  (Dc)  was  set  to  Dc  =  2rc  =  1.06  mm, 
corresponding  to  a  precipitation  rate  of  roughly  4  mm/h,  representing  a  light 
to  moderate  snow  rate  (e.g.,  Pruppacher  and  Klett  1980).  Because  snow 
particle  size  distributions  may  vary  greatly  for  a  given  precipitation  or  snow 
accumulation  rate,  it  is  difficult  to  make  generalizations  about  the  scattering 
characteristics  of  a  “typical”  snow.  Also,  due  to  the  use  of  the  AGAUS  Mie 
scattering  code,  snow  scattering  properties  are  being  modeled  using  a  spherical 
snowflake  assumption.  Winchester  et  al.  (1983)  note  “Since  the  experimental 
studies  have  shown  that  the  phase  functions  of  snow  crystals,  with  the  possible 
exception  of  graupel,  cannot  be  approximated  using  Mie  theory  computations 
for  spheres  with  either  equivalent  area  or  equivalent  volume.”  Thus,  admittedly 
this  is  a  poor  substitute  for  an  accurate  snow  model,  but  theoretical  techniques 
characterizing  non-spherical  scatterers  are  currently  inadequate  to  model  snow 
at  visible  and  infrared  wavelengths  on  non-super  computers.  Yet  extensions 
of  measurements  to  arbitrary  wavelength  are  impractical  without  some  model. 
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Compounding  these  difficulties  is  the  nature  of  the  behavior  of  the  refractive 
indices  of  ice.  These  are  temperature  and  humidity,  as  well  as  wavelength, 
dependent.  In  conclusion,  the  phase  function  presented  for  snow  can  at  best 
be  considered  an  approximation  to  the  actual  phase  function  for  snow,  and  at 
worst  a  simple  placeholder  for  a  future,  more  robust,  representation  of  variable 
snow  scattering  characterization. 

The  extinction  coefficients  for  other  precipitation  rates  may  be  calculated  using 
the  algorithm  given  in  XSCALE  (Fiegel  1994).  The  extinction  coefficients  for  the 
rain  and  snow  distributions  axe  presented  in  table  10. 


Figure  1.  Measured/modeled  snow  distribution 


2.5  Dust  Aerosol  Models 

Soil-derived  aerosols  are  an  important  component  of  the  total  atmospheric 
aerosol  content  in  certain  geographic  locations.  Reported  results  of  size 
distribution  measurements  for  these  aerosols  vary  widely.  However,  the  general 
consensus  is  that  the  dust  aerosols  follow  a  bimodal  lognormal  distribution. 
Empirical  data  (Patterson  and  Gillette  1977)  fits  this  type  of  distribution 
well,  and  dust  aerosols  may  be  produced  bj^  a  pulverization  process  in  the 
soil.  Epstein  (1947)  has  shown  that  such  processes  result  in  lognormal 
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Table  10.  Extinction  coefficients  (km  versus  wavelength  (//m)  for  various 
natural  aerosol  models 


Wavelength 

Rain 

Snow 

Dust  1 

(nm) 

Light 

Moderate 

Heavy 

0.55 

0.3664 

1.009 

1.561 

0.1179 

0.1018 

4.319 

1.06 

0.3673 

1.011 

1.564 

0.1180 

0.05754 

5.195 

3.00 

0.3697 

1.017 

1.572 

0.1184 

0.02799 

2.605 

3.50 

0.3705 

1.018 

1.574 

0.1186 

0.02557 

2.001 

4.00 

0.3711 

1.020 

1.576 

0.1186 

0.01907 

1.576 

4.50 

0.3716 

1.021 

1.578 

0.1187 

0.01513 

1.163 

5.00 

0.3721 

1.022 

1.579 

0.1188 

0.01276 

0.9553 

8.00 

0.3747 

1.028 

1.588 

0.1192 

0.00675 

0.6860 

8.50 

0.3751 

1.029 

1.589 

0.1193 

0.02951 

4.695 

9.00 

0.3755 

1.029 

1.590 

0.1193 

0.04410 

3.732 

9.50 

0.3758 

1.030 

1.591 

0.1194 

0.02062 

2.592 

10.00 

0.3761 

1.031 

1.592 

0.1194 

0.01598 

1.357 

10.50 

0.3760 

1.031 

1.592 

0.1195 

0.01235 

0.9214 

11.00 

0.3754 

1.030 

1.590 

0.1194 

0.009668 

0.7367 

11.50 

0.3745 

1.028 

1.588 

0.1194 

0.007398 

0.6250 

12.00 

0.3743 

1.028 

1.588 

0.1194 

0.004920 

0.5651 

Number 

density 

0.001951 

0.002736 

0.003165 

0.000040 

2480 

1258 

(cm“®) 

distributions.  The  bimodal  distribution  also  provides  a  better  fit,  as  empirical 
dust  distributions  appear  to  be  characterized  by  more  than  one  mode. 
Generally  (Patterson  and  Gillette  1977),  the  accumulation  or  small  mode 
appears  to  be  a  characteristic  of  dust  aerosols  under  all  conditions,  while  the 
coarse  or  large  mode  is  more  a  function  of  the  parent  soil  size  distribution.  The 
latter  component  usually  appears  only  under  conditions  of  moderate  to  heavy 
aerosol  dust  loading. 

The  parameters  for  light  and  heavy  aerosol  loading  (table  11)  were  taken  at 
various  locations,  predominantly  in  the  southwestern  United  States.  Analyses 
of  the  small  particle  mode  showed  (Jennings  et  al.  1978)  that  the  constituents 
were  primarily  ammonium  sulfate,  carbon,  calcite,  sodium  nitrate,  quartz,  and 
montmorillonite  for  both  distributions.  The  particles  contained  in  the  large 
mode  were  seen  to  settle  quickly,  in  both  light  and  moderate  cases,  as  the 
windspeed  diminished.  The  accumulation  mode  was  considered  to  be  comprised 
of  80  percent  quartz  and  20  percent  montmorillonite  by  mass  for  the  heavy 
aerosol  loading  and  80  percent  ammonium  sulfate  and  20  percent  carbon  by 
mass  for  the  light  aerosol  loading  case,  to  allow  the  distribution  to  be  more 
representative  of  varying  geographic  locations  (Gillespie  and  Lindberg  1992). 
The  resulting  distributions  are  representative  of  windblown  dust,  not  vehicular 
or  HE  debris. 
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Table  11.  Values  of  lognormal  particle  size  distribution  parameters  (mode  radius 
rg  and  geometric  standard  deviation  ag) 


Dust  Type 

Light  Loading 

Heavy  Loading  | 

Mode 

Small 

Large 

Small 

Large 

Species 

Ammonium  Sulfate 

Quartz 

Montmorillonite 

Quartz 

Bulk  Density 

1.769 

2.32 

2.5 

2.32 

(g/cm^) 

Number  Density 

1988 

3.79 

39.62 

0.1128 

(cm-3) 

Mass  loading 

16 

40 

1000 

10000 

Tg  (^m) 

0.05 

0.5 

0.5 

15 

2.0 

2.0 

2.25 

1.6 

Species 

Carbon 

Quartz  | 

Bulk  Density 

1.8 

2.32 

(g/cm^) 

Number  Density 

488.5 

1218.6 

(cm-3) 

Mass  loading 

4 

4000 

(^g/m®) 

rg  {^lm) 

0.05 

0.5 

£2 _ 

2.0 

1.6 

The  refractive  indices  for  quartz  were  interpolated  from  Weinman  and 
Peterson  (1969)  for  0.55  and  1.06  /zm,  from  Jennings  and  Gillespie  (1978)  for 
3.0  to  5.0  ^m,  and  from  Spitzer  and  Kleinman  (1961)  for  8.0  to  12.0/zm.  The 
refractive  indices  for  ammonium  sulfate  were  interpolated  from  the  work  of  Toon 
et  al.  (1976);  the  refractive  indices  of  carbon  were  interpolated  from  the  work  of 
Gillespie  and  Goedecke  (1989). 

The  refractive  indices  of  montmorillonite  at  0.55  and  1.06  /xm  were  interpolated 
from  the  work  of  Egan  and  Hilgeman  (1979);  for  wavelengths  greater  than  4.5 
/zm,  values  were  interpolated  from  the  work  of  Toon  et  al.  (1977). 

The  errors  introduced  by  the  interpolation  are  probably  small  because  of  the 
close  proximity  of  the  wavelengths  used  here  and  tabulated  in  the  aforementioned 
references.  However,  for  the  3.0-  to  4.5-^m  band,  the  refractive  indices  had 
to  be  interpolated  between  2.6  and  5.0  /um.  Because  montmorillonite  is  a 
clay  material  with  water  chemically  bonded  in  its  crystal  lattice  structure,  the 
refractive  indices  and  derived  quantities  should  be  used  with  extreme  caution  in 
this  wavelength  band  (3.0  to  4.5  /zm).  Tables  12  and  13  present  the  refractive 
indices  used  for  the  various  constituents. 

The  heavy  loading  dust  type  reflects  very  large  mode  radius  constituents 
associated  with  high  wind  speeds.  The  light  loading  would  be  the  case  normally 


Table  12.  Real  (n)  and  imaginary  (k)  indices  of  refraction  for  light  dust 
constituents  ammonium  sulfate  and  carbon 


Wavelength 

(//m) 

Ammonium  Sulfate 
n  k 

Carbon 
n  k 

0.55 

1.53 

1.0x10-^ 

2.0 

1.0 

1.06 

1.51 

2.1x10“® 

2.0 

1.0 

3.00 

1.36 

8.9x10-2 

2.2 

1.2 

3.50 

1.62 

1.4x10“^ 

2.2 

1.2 

4.00 

1.55 

1.7x10-2 

2.2 

1.2 

4.50 

1.50 

7.9x10-2 

2.2 

1.2 

5.00 

1.46 

7.0x10-2 

2.2 

1.2 

8.00 

1.31 

8.0x10-2 

3.0 

1.6 

8.50 

0.90 

2.7x10“^ 

3.0 

1.6 

9.00 

0.99 

1.7x10° 

3.0 

1.6 

9.50 

2.70 

6.1x10-2 

3.0 

1.6 

10.00 

2.19 

1.3x10-2 

3.0 

1.6 

10.50 

1.99 

6.0x10-2 

3.0 

1.6 

11.00 

1.90 

4.3x10-2 

3.0 

1.6 

11.50 

1.83 

2.8x10-2 

3.0 

1.6 

12.00 

1.80 

2.0x10-2 

3.0 

1.6 

Table  13.  Real  (n)  and  imaginary  (k)  indices  of  refraction  for  heavy  dust 
constituents  quartz  and  montmorillonite.  Quartz  results  include  ordinary  and 
extraordinary  indices. 


■■ 

Quartz 

Montmorillonite  1 

Bl 

n 

k 

n 

k 

Ordnry 

Extrord 

Ordnry 

Extrord 

0.55 

1.546 

1.555 

10“^ 

10“^ 

1.524 

0.000673 

1.06 

1.534 

1.543 

10-^ 

10-" 

1.519 

0.00057 

3.00 

1.500 

1.500 

10“® 

10“® 

1.483 

0.00317 

3.50 

1.485 

1.485 

10-® 

10"® 

1.463 

0.00350 

4.00 

1.472 

1.476 

0.00013 

0.00014 

1.442 

0.00383 

4.50 

1.426 

1.432 

0.00066 

0.00073 

1.421 

0.00417 

5.00 

1.412 

1.419 

0.00079 

0.00091 

1.400 

0.00450 

8.00 

0.42984 

0.39076 

0.13829 

0.14379 

1.035 

0.125 

8.50 

0.11260 

0.08548 

1.25062 

1.21601 

0.754 

0.427 

9.00 

0.17463 

0.22905 

2.59701 

3.04158 

0.923 

0.869 

9.50 

4.51517 

3.90448 

0.39770 

0.23041 

1.750 

1.860 

10.00 

2.66527 

2.57228 

0.05190 

0.04402 

2.590 

0.625 

10.50 

2.23766 

2.20003 

0.02452 

0.02220 

1,970 

0.185 

11.00 

2.01345 

2.00007 

0.01736 

0.01588 

1.845 

0.245 

11.50 

1.83358 

1.84954 

0.01875 

0.01529 

1.765 

0.160 

12.00 

1.56521 

1.68256 

0.04694 

0.02369 

1.693 

0.128 
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considered.  Mixing  the  heavy  and  light  cases  should  simulate  intermediate 
condition  dust  cases. 

Table  13  provides  information  on  both  the  ordinary  and  extraordinary  indices 
of  refraction  for  quartz.  Because  quartz  is  an  optically  positive  uniaxial 
crystal  (Born  and  Wolf,  1975)  the  scattering  problem  is  divided  into  two  parts. 
Two-thirds  of  the  scattering  material  is  treated  using  the  ordinary  indices  of 
refraction.  The  remaining  third  of  the  material  is  treated  using  the  extraordinary 
indices. 

The  HE  dust  model  was  generated  using  the  empirical  results  of  field 
tests  (Pinnick  et  al.  1983)  taken  at  Huntsville,  Alabama;  and  Orogrande,  New 
Mexico.  The  results  were  empirically  fitted  to  a  bimodal  lognormal  curve 
with  the  following  parameters:  for  the  small  mode,  particle  concentration 
C  =  15930  number  density  N  =  200cm“^,  geometric  mean  radius 

Tg  =  0.5  fim,  geometric  standard  deviation  ag  —  2.6,  and  particle  bulk  density 
p  =  2.5g/cm^;  for  the  large  mode,  C  =  48680  /ig/m^,  N  =  O.OTcm"^, 
rg  =  22.5  /zm,  Gg  =  1.87,  and  p  =  2.5g/cm^.  The  refractive  indices  were 
talcen  from  the  work  of  Ivlev  and  Popova  (1973).  The  refractive  indices  are  a 
synthetic  spectra  chosen  because  no  consistent  set  of  measurements  covers  the 
wavelength  range  in  PFNDAT.  A  comparison  of  the  synthesized  spectrum  with 
the  measurements  from  Jennings  et  al.  (1978)  at  the  wavelengths  available  shows 
agreement.  Table  14  shows  that  the  values  taken  fall  between  the  minimum  and 
maximum  values  found  in  that  reference.  Table  15  lists  the  extinction  coefficients 
for  all  smoke  types  along  with  the  results  for  the  HE  dust  type  as  functions  of 
wavelength,  as  determined  from  runs  of  the  AGAUS  model. 

Table  14.  Derived  real  (n)  and  imaginary  {k)  indices  of  refraction  for  HE  dust 


A 

n 

k 

0.55 

1.65 

0.005 

1.06 

1.647 

0.0051 

3.00 

1.646 

0.076 

3.50 

1.655 

0.020 

4.00 

1.637 

0.018 

4.50 

1.620 

0.018 

5.00 

1.592 

0.018 

8.00 

1.269 

0.178 

8.50 

1.186 

0.600 

9.00 

1.650 

1.240 

9.50 

2.342 

0.600 

10.00 

2.140 

0.126 

10.50 

1.904 

0.078 

11.00 

1.751 

0.143 

11.50 

1.784 

0.331 

12.00 

1.756 

0.230 

Table  15.  Extinction  coefficients  (km  versus  wavelength  (//m)  for  various 
manmade  aerosol  models 


Wavelength 

HE  Dust 

White  Phosphorus 

HC 

(pm) 

17%  Rh 

50%  Rh 

90%  Rh 

0.55 

2.668 

4191. 

4282. 

3957. 

5367. 

3227. 

1.06 

2.806 

1708. 

1963. 

2329. 

3737. 

2601. 

3.00 

309.1 

449.6 

596.9 

1141. 

3.50 

414.8 

421.5 

342.1 

641.4 

383.6 

4.00 

287.2 

297.7 

187.3 

4.50 

255.3 

258.5 

208.2 

143.8 

5.00 

2.918 

148.9 

8.00 

1.944 

548.2 

519.6 

38.61 

8.50 

2.399 

421.3 

421.7 

31.81 

53.5 

9.00 

25.31 

49.6 

9.50 

2.944 

516.2 

248.4 

22.42 

49.5 

10.00 

2.892 

236.2 

19.34 

57.4 

10.50 

2.713 

361.5 

17.27 

67.7 

11.00 

248.9 

226.5 

84.5 

11.50 

2.515 

136.3 

157.1 

14.92 

12.00 

2.432 

117.4 

138.7 

182.6 

12.49 

141.8 

Number 

density 
(10®  cm-3) 

0.0002 

5.667 

4.566 

2.239 

8.261 

1.399 

2.6  Aerosol  Smoke  Models 

The  phase  functions  for  inventory  smokes  can  be  calculated  almost  exactly 
because  the  particles  are  nearly  spherical.  Discrepancies  between  theory  and 
measurement  can  be  attributed  to  uncertainties  in  the  particle  size  spectrum 
or  complex  refractive  indices.  Experiments  (Jennings  and  Gillespie  1978)  have 
shown  that  the  particle  size  spectrum  is  closely  approximated  by  a  lognormal 
distribution.  Reference  to  the  mass  loading  or  mass  concentration  (C)  of  the 
particulate  material  (equation  4)  rather  than  the  number  density  is  conventional 
in  smoke  applications.  Since  the  bulk  density  (p)  for  water  is  1  g/cm^,  p  for  a 
smoke  is  also  numerically  equal  to  the  particulate  specific  gravity. 

Table  16  lists  the  parameters  considered  representative  of  inventory  smokes  and 
includes  the  mass  median  diameter  (MMD)  often  used  in  the  literature  in  place 
of  Vg.  The  two  are  related  as 

In(MMD)  =  In  2rg  +  3  In^  ag ,  (22) 

where  Vg  and  ag  are  listed  in  table  16. 

Table  16  shows  that  mass  concentration  C  was  arbitrarily  set  to  10®  pg/m^. 
The  magnitude  of  C  has  no  effect  on  the  phase  function  or  the  mass  extinction 
coefficient  and  only  linearly  scales  the  volume  extinction  coefficient  /Jj. 
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Table  16.  Representative  parameters  for  determining  phase  functions  of 
inventory  smokes  at  various  relative  humidities 


Aerosol  Species: 

Relative  Humidity  (%): 

Whii 

te  Phosphorus 

Fog  Oil 

HC 

17%  Rh 

50%  Rh 

90%  Rh 

Geometric  mean  (^m),  rg 
Width  parameter,  ag 

Bulk  density  (g/cm^),  p 

Mass  loading  (//g/m^),  C 

Mass  median  diameter,  MMD 

0.241 

1.450 

1.617 

10® 

0.729 

0.269 

1.450 

1.443 

10® 

0.814 

0.365 

1.450 

1.178 

10® 

1.104 

0.190 

1.800 

0.890 

10® 

0.575 

0.422 

1.450 

1.220 

10® 

1.338 

Table  16  lists  the  fog  oil  particle  spectrum  parameters  appropriate  for  fog 
oil  dissemination  by  current  military  generators  designed  to  produce  particles 
most  efficient  for  obscuration  at  the  visible  wavelengths  (Carlon  et  al.  1977). 
Other  experimental  generators  may  produce  larger  particles.  WP  and  HC  have 
parameters  listed  for  specific  values  of  relative  humidity.  Hygroscopic  growth  has 
been  modeled  for  these  conditions  by  semiempirical  relations  (Frickel  et  al.  1979; 
Rubel  1978).  Other  evidence  (Farmer  1980)  shows  that  at  high  humidities 
(greater  than  75  percent  relative  humidity)  a  bimodal  particle  size  spectrum 
may  be  expected  that  would  be  most  pronounced  for  WP  smoke. 

Tabulated  real  and  imaginary  refractive  indices  have  also  been  provided. 
Reliable  experimental  measurements  would  be  preferred,  but  such  measurements 
axe  usually  impossible  (Weast  and  Astle  1980)  because  of  the  complex 
reaction  products  formed  in  producing  smoke.  Table  17  lists  the  utilized 
values  for  selected  wavelengths  from  the  visible  through  the  infrared. 
Weast  and  Astle  (1980)  is  the  primary  reference.  They  derived  coefficients  based 
on  laboratory  measurements  performed  on  the  major  constituents  —  phosphoric 
acid,  H3PO4,  and  zinc  chloride,  ZnCl2,  for  WP  and  HC  scatterers,  respectively, 
at  various  humidity  levels. 

Fog  oil  smoke  is  not  considered  hygroscopic,  so  only  a  single  data  set  is  used. 
In  the  visible,  the  imaginary  index  for  fog  oil  is  so  small  as  to  be  beyond 
instrumental  sensitivity;  it  can  be  considered  negligible  for  most  applications. 
This  small  value  for  k  leads  to  a  single-scattering  albedo  of  nearly  unity,  implying 
that  extinction  is  entirely  due  to  scattering. 

Table  18  compares  the  average  mass  extinction  coefficients  as  computed 
by  AGAUS  versus  laboratory  experimental  results  (Weast  and  Astle  1980; 
Frickel  et  al.  1979;  Rubel  1978;  Farmer  1980)  for  several  spectral  bands  of 
interest.  Since  the  measured  results  represent  band  averages,  a  typical  cloud 
thickness  was  assumed  (R  =  0.01  km)  and  results  were  computed  by  averaging 
the  computed  transmission  through  1  g/m^  density  aerosols  via. 


A:  =  —  In 


Wi  exp(— fcjJZ)  ' 


(23) 
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Table  17.  Real  (n)  and  imaginary  (k)  indices  of  refraction  for  the  smoke  aerosol 
models  at  indicated  relative  humidities 


H3PO4  (WP)  17%  RH  H3PO4  (WP)  50%  RH  H3PO4  (WP)  90%  RH 


aSmaM 

n 

k 

n 

k 

n 

k 

0.55 

1.438 

0.001 

1.412 

0.0008 

1.357 

0.0003 

1.06 

1.414 

0.008 

1.399 

0.0057 

1.348 

0.0018 

3.00 

1.278 

0.104 

1.301 

0.133 

1.350 

0.2290 

3.50 

1.356 

0.178 

1.363 

0.150 

1.389 

0.0522 

4.00 

1.338 

0.141 

1.382 

0.118 

1.360 

0.0393 

4.50 

1.417 

0.150 

1.403 

0.127 

1.354 

0.0481 

5.00 

1.399 

0.119 

1.387 

0.101 

1.344 

0.0395 

8.00 

1.287 

0.622 

1.288 

0.524 

1.290 

0.184 

8.50 

1.421 

0.557 

1.383 

0.480 

1.310 

0.172 

9.00 

1.396 

0.615 

1.374 

0.519 

1.296 

0.186 

9.50 

1.462 

0.807 

1.510 

0.665 

1.304 

0.242 

10.00 

1.720 

0.827 

1.636 

0.697 

1.346 

0.248 

10.50 

1.793 

0.826 

1.691 

0.699 

1.340 

0.259 

11.00 

2.125 

0.768 

1.962 

0.671 

1.400 

0.272 

11.50 

2.080 

0.404 

1.920 

0.360 

1.368 

0.209 

12.00 

1.951 

0.329 

1.810 

0.307 

1.324 

0.232 

RH 

Fog  Oil,  50%  RH 
n  k 

ZnCl2  (HC),  85%  RH 
n  k 

0.55 

1.475 

0.000002 

0.000 

1.06 

1.474 

0.000006 

BR:« 

0.000 

3.00 

1.466 

0.000337 

B|E!^R 

0.227 

3.50 

1.518 

0.0466 

0.021 

4.00 

1.482 

0.000701 

1.405 

0.005 

4.50 

1.479 

0.000504 

1.382 

0.016 

5.00 

1.476 

0.000357 

1.376 

0.018 

8.00 

1.485 

0.00491 

1.348 

0.037 

8.50 

1.480 

0.00514 

1.336 

0.040 

9.00 

1.480 

0.00407 

1.321 

0.041 

9.50 

1.478 

0.00504 

1.300 

0.045 

10.00 

1.479 

0.00509 

1.279 

0.057 

10.50 

1.479 

0.00557 

1.253 

0.072 

11.00 

1.479 

0.00467 

1.229 

0.095 

11.50 

1.479 

0.00714 

1.204 

0.128 

12.00 

1.478 

0.00620 

1.186 

0.174 

where  the  Wi  represent  weight  factors  for  each  spectral  band.  A  simple  approach 
sets  the  weight  factors  in  the  first  and  last  bands  to  1/2,  and  the  remaining 
weights  to  1.  This  was  the  method  used  to  produce  the  model  results  of  table  18. 
This  table  shows  that  all  the  comparisons  are  reasonable.  Disparities  are  no 
larger  than  those  found  among  various  experiments  throughout  the  above  cited 
literature.  Comparison  of  results  for  WP  in  the  8  to  12- /xm  region  are  sometimes 
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taken  as  evidence  (Milham  et  al.  1977)  that  secondary  reaction  products  are 
significant  for  WP  smokes.  Results  for  fog  oil  at  visible  wavelengths  may 
be  due  to  the  use  of  a  single  wavelength  for  the  modeled  results.  Results 
in  the  longwave  IR  band  may  reflect  different  assumptions  regarding  sources. 
Equation  (23)  assumes  a  flat  source  spectrum.  The  measurements  would  rely 
on  the  temperature  of  the  background  medium.  It  is  probably  significant 
that  the  experimental  data  were  obtained  by  the  vapor  condensation  method 
rather  than  by  pyrotechnic  dissemination.  Disparities  were  noted  before  by 
Pinnick  and  Jennings  (1980). 

Table  18.  Comparison  of  theoretical  (from  the  AGAUS  model)  and 
experimentally  measured  (Expt)  mass  extinction  coefficients  (m^/g)  at  50 
percent  relative  humidity  for  various  smoke  aerosols 


A  (^m) 

WP 

Model 

Expt 

Fog  Oil 

Model  Expt 

HC 

Model 

Exp 

visible* 

4.282 

3.940 

5.367 

7.730 

3.227 

4.579 

1.06 

1.963 

1.410 

3.737 

3.500 

2.601 

2.040 

3-5 

0.284 

0.290 

0.262 

0.270 

0.193 

0.190 

8-12 

0.284 

0.366 

0.021 

0.014 

0.068 

0.052 

*0.55  jJLVCi  for  model;  0.4  to  0.7  ^m  for  experiment 


Of  further  interest  to  the  usage  of  phase  functions  for  smoke  aerosols  are  the 
single  scattering  albedos  (w)  averaged  over  various  wavebands.  Table  19  lists  the 
average  single-scattering  albedos  for  the  inventory  smokes  in  four  spectral  regions 
of  interest.  Due  to  the  usage  of  updated  fog  oil  imaginary  refractive  indices 
approximately  an  order  of  magnitude  lower  than  those  used  in  the  EOSAEL87 
version  of  PFNDAT,  the  fog  oil  single  scattering  albedos  are  considerably  higher 
in  the  IR  bands  than  previously  reported.  The  significancy  of  this  updated 
finding  is  that  scattering  becomes  more  significant  even  at  IR  wavelengths. 

Table  19.  Average  single-scattering  albedo  for  the  inventory  smokes  as  calculated 
by  AGAUS 


A  (//m) 

WP 

Fog  Oil 

HC 

0.55 

0.995 

>  0.999 

>  0.999 

1.06 

0.964 

>  0.999 

>  0.999 

3-5 

0.155 

0.916 

0.745 

8-12 

0.017 

0.652 

0.081 

The  phase  function  does  not  depend  on  number  density,  but  the  volume 
extinction  coefficient  does.  Thus  table  16  lists  the  parameter  values  used 
to  compute  the  extinction  coefficients  for  the  various  smokes.  Because  the 


smoke  density  is  a  definite  function  of  time,  a  method  for  reassessing  the 
extinction  coefficient  at  different  times  is  necessary.  This  reassessment  may 
be  accomplished  as  follows:  In  the  notation  used  here,  the  transmission  T  is. 


T  =  exp  ^ 

(24) 

=  exp“"'*^^ 

(25) 

O 

a 

II 

(26) 

=  Ns, 

(27) 

where 

•  L  is  the  path  length  (cm), 

•  13 f  is  the  volume  extinction  coefficient  (cm“^), 

•  Of  is  the  mass  extinction  coefficient  (cm^/g),  and 

•  s  is  the  extinction  cross  section  per  particle  (cm^). 

Thus  the  volume  extinction  coefficient  13^  can  be  scaled  as  a  function  of  time  if 
N  (or  the  quantity  OeC)  is  known  as  a  function  of  time.  In  addition,  L  will  vary 
with  time  according  to  statistical  variations,  elapsed  time  since  the  smoke  event 
began,  and  the  observer’s  geometry  with  respect  to  the  cloud.  The  COMBIC 
module  contained  in  EOSAEL  provides  mean  estimates  for  the  quantities  T  and 
L  for  specified  LOS  through  various  geometries  of  smoke  clouds. 

2.7  EOSAEL92  Improvements 

The  distinctions  between  the  current  version  of  PEN  DAT  (part  of  the  1992  release 
of  EOSAEL  (EOSAEL92))  and  the  1987  release  of  PFNDAT  (part  of  EOSAEL87, 
the  EOSAELversion  released  in  1987)  are  significant.  This  version  has  expanded 
the  visible  band  calculations  to  include  results  at  0.05-/Ltm  resolution  from  0.35 
through  0.75  fj.m.  This  version  of  the  documentation  also  contains  sufficient 
information  for  the  user  to  duplicate  the  calculations  made  using  the  AGAUS 
model  via  a  different  Mie  scattering  routine.  The  version  of  AGAUS  utilized 
was  updated  to  include  a  continuing  fraction  expansion  technique  proposed  by 
Lentz  (1976).  This  approach  allows  for  a  more  precise  computation  of  forward 
scattering  effects.  The  forward  peaks  of  some  of  the  larger  aerosols  (most  notably 
rain  and  snow)  have  increased  significantly,  better  reflecting  the  true  forward 
scattering  effects. 

We  have  updated  the  index  of  smoke-aerosol  refraction  data  used  for  WP  and 
fog  oil  to  reflect  more  recent  data  (Hoock  and  Sutherland  1993).  The  reported 
infrared  imaginary  indices  of  refraction  for  fog  oil  were  adjusted  (following 
remeasurement)  from  the  values  used  in  the  original  PFNDAT.  Minor  errors 
in  the  text  were  also  corrected,  as  well  as  inconsistencies  in  notation  used  in  the 
original  document. 
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3.  USER’S  GUIDE 

3.1  Introduction 

The  aerosol  phase  function  data  files  PFNDAT .  nnn  (where  nnn  varies  from  1  to  57) 
are  accessed  by  modules  ASCAT  and  FCLOUD  for  a  given  aerosol  distribution 
and  wavelength.  The  selected  phase  function  is  renormalized  by  subroutine 
PFUNC  for  the  MSCAT  module  (an  EOSAEL  Multiple  SCATtering  routine)  and 
by  subroutine  PFN  for  the  FCLOUD  module  so  its  normalization  is  compatible 
with  the  calling  module.  ASCAT  requires  no  renormalization. 

Computer  code  AGAU5  and  a  geometrical  optics  version  of  code  AGSNOW  were 
used  to  construct  the  aerosol  phase  function  data  base.  Table  20  lists  the  38 
different  distributions  contained  in  the  database. 

3.2  Usage 

To  use  PFNDAT  in  one  of  the  aforementioned  codes  requires  the  assignment  of 
the  FORTRAN  unit  number,  IPHFUN,  within  the  codes.  This  assignment  depends 
on  the  specific  location  of  the  database  within  the  user’s  computer  system. 
The  PFNDAT  files  are  provided  in  two  primary  forms.  In  one,  the  individual 
scattering  types  are  located  in  separate  output  files.  In  the  other,  the  contents 
of  the  individual  files  are  concatenated  to  produce  a  large  master  file. 

If  users  wish  to  construct  their  own  phase  function(s)  for  use  by  one  of  the 
programs,  this  file  must  be  formatted  as  explained  in  the  structure  section  below. 
Program  AGAUS,  supplied  as  an  ancillary  code  to  EOSAEL,  may  be  used  to 
automatically  construct  a  file  compatible  with  EOSAEL  usage.  Instructions  for 
the  use  of  AGAUS  are  found  as  comments  at  the  beginning  of  that  code. 

3.3  Structure 

The  phase  function  database  comprises  a  series  of  ASCII  files,  one  for  each  of 
the  aerosol  identifiers  listed  in  table  20.  The  files  are  called  PFNDAT .  001  through 
PFNDAT. 057.  Each  identifier  nnn  is  associated  with  a  file  PFNDAT. nnn.  Each  file 
begins  with  65  discrete  angles  between  0°  and  180°;  the  number  of  angles  is 
the  current  dimension  size  (65)  of  the  pertinent  arrays  in  the  EOSAEL  module 
cited  above.  The  remainder  of  the  file  contains  sets  of  phase  function  results  at 
each  wavelength.  Each  set  is  composed  of  a  one-line  preamble  followed  by  the 
angular  phase  function  data.  The  preamble  is  a  record  containing  the  number 
of  angular  data  items  (NANG,  65  in  all  cases),  a  phase  function  identifier 
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Table  20.  Phase  function  data  base  for  EOSAEL92 


Index 

Distribution  Type 

%  Rel.  Hmdty. 

1. 

Maritime 

0 

2. 

Maritime 

50 

3. 

Maritime 

70 

4. 

Maritime 

80 

5. 

Maritime 

90 

6. 

Maritime 

95 

7. 

Maritime 

98 

8. 

Maritime 

99 

9. 

Urban 

0 

10. 

Urban 

50 

11. 

Urban 

70 

12. 

Urban 

80 

13. 

Urban 

90 

14. 

Urban 

95 

15. 

Urban 

98 

16. 

Urban 

99 

17. 

Rural 

0 

18. 

Rural 

50 

19. 

Rural 

70 

20. 

Rural 

80 

21. 

Rural 

90 

22. 

Rural 

95 

23. 

Rural 

98 

24. 

Rural 

99 

25. 

Fog  (heavy  advection) 

NA 

26. 

Fog  (moderate  radiation) 

NA 

27. 

Rain  (drizzle) 

NA 

28. 

Rain  (widespread) 

NA 

29. 

Rain  (thunderstorm) 

NA 

30. 

Snow 

NA 

31.-49. 

(Reserved  for  future  use) 

50. 

Dust  (light  loading) 

NA 

51. 

Dust  (heavy  loading) 

NA 

52. 

High  explosive  (HE)  dust 

NA 

53. 

WP  smoke 

17 

54. 

WP  smoke 

50 

55. 

WP  smoke 

90 

56. 

Fog  oil 

50 

57. 

HC  smoke 

85 

(0  =  user  supplied),  wavelength  of  this  set  (micrometers),  the  albedo  for  single 
scattering,  and  the  extinction  and  scattering  coefficients  in  inverse  kilometers. 
Subsequent  to  the  preamble  are  values  of  the  phase  function  at  each  angle.  For 
most  scattering  species  there  will  be  32  such  sets  of  angular  data.  For  the  smokes. 


only  20  wavelength  sets  are  provided,  since  data  was  unavailable  at  wavelengths 
beyond  12  yun..  Table  21  is  a  schematic  example  of  the  structure  for  PFNDAT. 

The  PFNDAT  phase  function  data  is  formatted  such  that  there  is  always  one 
more  data  item  than  the  number  of  angular  results  called  for.  The  la.st  value  is 
always  set  to  999.99,  which  is  used  by  the  PFUNC  routine  to  determine  the  end 
of  the  angular  data.  If  a  user- specified  scattering  species  is  used  that  has  fewer 
than  65  angles,  then  a  single  additional  value  must  be  included  in  the  file  to 
indicate  the  end  of  each  set  of  angular  varying  phase  function  data. 

Subroutine  PFUNC  counts  the  number  of  angles,  looking  for  a  value  of  999.99,  and 
will  compare  this  internally  counted  number  with  the  value  of  N AN G,  the  total 
number  of  input  angles.  Should  the  two  numbers  disagree,  an  error  message 
is  printed  and  execution  is  halted.  The  fog,  rain,  and  snow  distributions  use 
a  different  set  of  angles  that  are  included  in  subroutine  PFUNC;  the  angles  are 
automatically  invoked  by  using  the  phase  function  identifier  as  a  switch. 


39 


Table  21.  Structure  for  Aerosol  Phase  Function  Data  File  PFNDAT.nnn 


^56 

^2 

^65 

Oil 

999.99 

NANG 

nnn 

/?ei 

0s 

P(6ii,  Ai,nnn) 

P(^2,  Ai,nnn) 

P(^6,  Ai,nnn) 

P(06i,  Ai,nnn) 

P(^2,  Ai,nnn) 

.  .  . 

-P(^65,  Ai,nnn) 

NANG 

nnn 

^2 

Pex 

0s 

P(6'i,  A2,nnn) 

P(^2,  A2,nnn) 

.  .  . 

P(^6,  A2,nnn) 

P(^6i,A2,nnn) 

P(^2,  A2,nnn) 

.  .  . 

P(^65,  A2,nnn) 

:  :  :  :  :  : 

NANG 

nnn 

^max 

^ex 

Co 

P(^i,  A^ai,nnn) 

■P(^2  ?  ^max-i  nnn) 

.  .  . 

-P(^6  9  ^max'f  nnn) 

•)  ^max'i  nnn) 

^max>i  nnn) 

.  .  . 

P(^65y  ^maxi  nnn) 

0i  = 

NANG  = 
nnn  = 
A  = 


/^ex  — 

P(^i,  A/,nnn)  = 


discrete  angles  (degrees) 

number  of  discrete  angles 

phase  function  identifier  from  table  20 

wavelength  (micrometers) 

albedo  for  single  scattering 

extinction  coefiicient  (km“^ ) 

scattering  coefficient  (km“^) 

the  value  of  the  phase  function  at  angle  i,  wavelength  i, 
and  aerosol  type  identifier  nnn 
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ACRONYMS  and  ABBREVIATIONS 


ASL 

AFGL 

ARL 

AGAUS 

AGSNOW  - 

AS  CAT 

EOSAEL  - 

EOSAEL92  - 

FCLOUD  - 

HC 

HE 

LOS 

MSCAT 

WP 


The  U.S.  Army  Atmospheric  Sciences  Laboratory 

The  Air  Force  Geophysics  Laboratory 

The  U.S.  Army  Research  Laboratory 

August  Miller’s  Mie  Scattering  Code 

A  special  version  of  AGAUS  for  treating  SNOW  cases 

An  Approximate  multiple  Scattering  model  within  EOSAEL 

The  Electro- Optical  Systems  Atmospheric  Effects  Library 

The  1992  Release  of  EOSAEL 

A  Finite  Cloud  Transmission  model  within  EOSAEL 
Hexachloroethane  smoke  munition 
High  Explosive  artillery  munition 
Line  of  Sight 

A  Multiple  Scattering  model  within  EOSAEL 
White  Phosphorus  smoke  munition 
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Appendix  A 

INDICES  OF  REFRACTION 


The  refractive  index  information  for  the  various  haze  aerosol  constituents  is 
included  in  this  appendix  for  consistency  and  completeness.  The  rural  aerosol 
is  composed  of  small  and  large  rural  aerosol  components.  The  urban  aerosol  is 
composed  of  small  and  large  urban  aerosol  components.  The  maritime  aerosol 
is  composed  of  the  small  rural  and  oceanic  aerosol  components.  Care  is  needed 
in  composing  the  correct  constituent  components  of  each  aerosol.  The  aerosol 
extinction  coefficient  information  listed  in  the  main  text  is  produced  by  using 
the  specific  particle  size  distribution  information  contained  in  the  tables  in  the 
main  text  and  the  index  of  refraction  information  contained  in  this  appendix  as 
input  to  the  AGAUS  model  in  a  Mie  scattering  calculation. 
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Table  A-1.  Index  of  refraction  as  a  function  of  relative  humidity  (0-80  %)  for 
small  urban  aerosols 


Table  A-2.  Index  of  refraction  as  a  function  of  relative  humidity  (90-99  %)  for 
small  urban  aerosols 


Relative  Humiditvf%)  i 

A 

90 

95 

98 

99  1 

im) . 

n 

k 

n 

k 

n 

k 

n 

k 

0.3371 

1.394 

0.0210 

1.375 

0.0131 

1.362 

0.00716 

1.356 

0.00480 

0.4000 

1.389 

0.0206 

1.370 

0.0128 

1.356 

0.00701 

1.350 

0.00471 

0.4880 

1.386 

0.0202 

1.367 

0.0125 

1.352 

0.00687 

1.347 

0.00461 

0.5145 

1.385 

0.0202 

1.366 

0.0125 

1.351 

0.00687 

1.346 

0.00461 

0.5500 

1.384 

0.0199 

1.365 

0.0124 

1.350 

0.00676 

1.345 

0.00454 

0.6328 

1.384 

0.0194 

1.364 

0.0121 

1.350 

0.00662 

1.344 

0.00444 

0.6943 

1.383 

0.0196 

1.363 

0.0122 

1.349 

0.00666 

1.343 

0.00447 

0.8600 

1.379 

0.0201 

1.360 

0.0125 

1.346 

0.00686 

1.341 

0.00461 

1.0600 

1.377 

0.0212 

1.358 

0.0132 

1.343 

0.00721 

1.338 

0.00484 

3.0000 

1.386 

0.240 

1.380 

0.252 

1.376 

0.261 

1.374 

0.265 

3.5000 

1.420 

0.0324 

1.413 

0.0237 

1.407 

0.0172 

1.405 

0.0147 

4.0000 

1.383 

0.0295 

1.371 

0.0201 

1.362 

0.0131 

1.358 

0.0103 

4.5000 

1.369 

0.0379 

1.355 

0.0286 

1.345 

0.0218 

1.341 

0.0190 

5.0000 

1.364 

0.0376 

1.349 

0.0281 

1.338 

0.0210 

1.334 

0.0182 

7.9000 

1.311 

0.0650 

1.304 

0.0533 

1.300 

0.0445 

1.298 

0.0410 

8.2000 

1.281 

0.0723 

1.283 

0.0582 

1.284 

0.0478 

1.285 

0.0436 

8.5000 

1.319 

0.0885 

1.303 

0.0689 

1.292 

0.0543 

1.287 

0.0485 

9.0000 

1.478 

1.12 

1.396 

0.0851 

1.335 

0.0646 

1.311 

0.0565 

9.5000 

1.392 

0.0925 

1.336 

0.0743 

1.294 

0.0607 

1.277 

0.0554 

10.0000 

1.359 

0.0896 

1.306 

0.0750 

1.266 

0.0640 

1.250 

0.0597 

10.5910 

1.315 

0.0986 

1.264 

0.0868 

1.225 

0.0780 

1.210 

0.0745 

11.0000 

1.290 

0.119 

1.238 

0.110 

1.200 

0.104 

1.184 

0.102 

11.5000 

1.262 

0.154 

1.211 

0.149 

1.172 

0.146 

1.157 

0.145 

12.5000 

1.251 

0.247 

1.203 

0.251 

1.167 

0.255 

1.152 

0.256 

14.0000 

1.314 

0.337 

1.275 

0.350 

1.246 

0.359 

1.234 

0.363 

15.0000 

1.348 

0.379 

1.319 

0.388 

1.297 

0.394 

1.288 

0.397 

18.0000 

1.542 

0.398 

1.497 

0.408 

1.464 

0.416 

1.450 

0.420 

20.0000 

1.608 

0.383 

1.560 

0.387 

1.524 

0.390 

1.509 

0.391 

25.0000 

1.637 

0.364 

1.597 

0.361 

1.567 

0.359 

1.555 

0.358 

30.0000 

1.639 

0.355 

1.606 

0.345 

1.581 

0.337 

1.571 

0.334 

35.0000 

1.644 

0.379 

1.602 

0.363 

1.570 

0.351 

1.558 

0.346 

40.0000 

1.639 

0.436 

1.594 

0.417 

1.560 

0.402 

1.546 

0.397 
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Table  A-3.  Index  of  refraction  as  a  function  of  relative  humidity  (0-80  %)  for 
large  urban  aerosols 


Relative  Humiditv(%)  1 

A 

0 

50 

70 

80  1 

(urn) 

n 

k 

n 

k 

n 

k 

n 

k 

0.3371 

1.574 

0.0987 

1.556 

0.0908 

1.479 

0.0580 

1.428 

0.0323 

0.4000 

1.574 

0.0967 

1.555 

0.0890 

1.477 

0.0568 

1.416 

0.0316 

0.4880 

1.574 

0.0947 

1.555 

0.0871 

1.475 

0.0556 

1.413 

0.0310 

0.5145 

1.574 

0.0947 

1.555 

0.0871 

1.475 

0.0556 

1.413 

0.0310 

0.5500 

1.574 

0.0933 

1.555 

0.0858 

1.474 

0.0548 

1.412 

0.0305 

0.6328 

1.574 

0.0913 

1.555 

0.0840 

1.474 

0.0536 

1.411 

0.0299 

0.6943 

1.574 

0.0918 

1.555 

0.0845 

1.474 

0.0539 

1.411 

0.0300 

0.8600 

1.566 

0.0946 

1.547 

0.0871 

1.468 

0.0556 

1.407 

0.0310 

1.0600 

1.566 

0.0994 

1.547 

0.0815 

1.467 

0.0584 

1.405 

0.0325 

3.0000 

1.442 

0.123 

1.436 

0.135 

1.412 

0.185 

1.394 

0.223 

3.5000 

1.495 

0.117 

1.488 

0.109 

1.456 

0.0728 

1.431 

0.0448 

4.0000 

1.501 

0.122 

1.489 

0.112 

1.439 

0.0733 

1.400 

0.0429 

4.5000 

1.508 

0.129 

1.494 

0.119 

1.435 

0.0811 

1.390 

0.0511 

5.0000 

1.506 

0.131 

1.492 

0.121 

1.431 

0.0818 

1.384 

0.0511 

7.9000 

1.372 

0.180 

1.366 

0.168 

1.340 

0.120 

1.320 

0.0817 

8.2000 

1.263 

0.210 

1.265 

0.196 

1.272 

0.134 

1.278 

0.0922 

8.5000 

1.470 

0.280 

1.455 

0.261 

1.391 

0.180 

1.341 

0.116 

9.0000 

2.276 

0.381 

2.194 

0.354 

1.857 

0.240 

1.594 

0.151 

9.5000 

1.945 

0.270 

1.889 

0.252 

1.655 

0.177 

1.473 

0.118 

10.0000 

1.881 

0.233 

1.828 

0.219 

1.607 

0.158 

1.435 

0.111 

10.5910 

1.818 

0.214 

1.767 

0.202 

1.554 

0.153 

1.388 

0.115 

11.0000 

1.798 

0.199 

1.746 

0.191 

1.532 

0.157 

1.364 

0.130 

11.5000 

1.765 

0.198 

1.714 

0.194 

1.501 

0.175 

1.335 

0.160 

12.5000 

1.724 

0.201 

1.676 

0.206 

1.476 

0.225 

1.320 

0.240 

14.0000 

1.700 

0.216 

1.661 

0.229 

1.498 

0.280 

1.370 

0.320 

15.0000 

1.638 

0.294 

1.609 

0.303 

1.486 

0.339 

1.390 

0.367 

18.0000 

1.982 

0.293 

1.938 

0.304 

1.751 

0.348 

1.606 

0.383 

20.0000 

2.080 

0.346 

2.032 

0.350 

1.832 

0.365 

1.676 

0.378 

25.0000 

2.028 

0.392 

1.988 

0.389 

1.823 

0.377 

1.693 

0.368 

30.0000 

1.965 

0.455 

1.932 

0.445 

1.794 

0.402 

1.687 

0.369 

35.0000 

2.057 

0.538 

2.015 

0.522 

1.840 

0.455 

1.704 

0.402 

40.0000 

2.084 

0.624 

2.038 

0.605 

1.850 

0.525 

1.704 

0.463 

Table  A-4.  Index  of  refraction  as  a  function  of  relative  humidity  (90-99  %)  for 
large  urban  aerosols 


Relative  Humiditv(%)  1 

A 

90 

95 

98 

99  1 

inm) 

n 

k 

n 

k 

n 

k 

n 

k 

0.3371 

1.387 

0.0179 

1.368 

0.00982 

1.354 

0.00395 

1.349 

0.00193 

0.4000 

1.382 

0.0176 

1.362 

0.00962 

1.348 

0.00387 

1.344 

0.00189 

0.4880 

1.378 

0.0172 

1.359 

0.00942 

1.345 

0.00379 

1.340 

0.00185 

0.5145 

1.378 

0.0172 

1.358 

0.00942 

1.344 

0.00379 

1.339 

0.00185 

0.5500 

1.377 

0.0170 

1.357 

0.00928 

1.343 

0.00374 

1.338 

0.00182 

0.6328 

1.376 

0.0166 

1.356 

0.00908 

1.342 

0.00366 

1.337 

0.00178 

0.6943 

1.375 

0.0167 

1.355 

0.00913 

1.341 

0.00368 

1.336 

0.00179 

0.8600 

1.372 

0.0172 

1.353 

0.00941 

1.338 

0.00379 

1.334 

0.00185 

1.0600 

1.370 

0.0181 

1.350 

0.00989 

1.336 

0.00399 

1.331 

0.00194 

3.0000 

1.384 

0.245 

1.378 

0.257 

1.374 

0.266 

1.372 

0.269 

3.5000 

1.417 

0.0290 

1.409 

0.0201 

1.404 

0.0137 

1.402 

0.0115 

4.0000 

1.378 

0.0259 

1.366 

0.0162 

1.357 

0.00929 

1.354 

0.00688 

4.5000 

1.364 

0.0344 

1.350 

0.0249 

1.339 

0.0180 

1.335 

0.0156 

5.0000 

1.358 

0.0339 

1.343 

0.0241 

1.332 

0.0171 

1.329 

0.0147 

7.9000 

1.308 

0.0605 

1.302 

0.0484 

1.297 

0.0398 

1.296 

0.0368 

8.2000 

1.282 

0.0669 

1.284 

0.0525 

1.285 

0.0421 

1.286 

0.0385 

8.5000 

1.313 

0.0809 

1.297 

0.0608 

1.286 

0.0464 

1.282 

0.0414 

9.0000 

1.446 

0.102 

1.363 

0.0738 

1.303 

0.0536 

1.282 

0.0466 

9.5000 

1.371 

0.0855 

1.313 

0.0668 

1.271 

0.0534 

1.257 

0.0488 

10.0000 

1.339 

0.0840 

1.264 

0.0689 

1.245 

0.0581 

1.231 

0.0544 

10.5910 

1.295 

0.0941 

1.243 

0.0820 

1.265 

0.0733 

1.191 

0.0703 

11.0000 

1.270 

0.115 

1.217 

0.107 

1.179 

0.101 

1.166 

0.0988 

11.5000 

1.242 

0.152 

1.190 

0.148 

1.152 

0.144 

1.138 

0.143 

12.5000 

1.232 

0.249 

1.133 

0.253 

1.147 

0.257 

1.135 

0.258 

14.0000 

1.299 

0.342 

1.259 

0.355 

1.230 

0.364 

1.220 

0.367 

15.0000 

1.337 

0.382 

1.307 

0.391 

1.285 

0.398 

1.277 

0.400 

18.0000 

1.525 

0.402 

1.475 

0.413 

1.445 

0.421 

1.434 

0.423 

20.0000 

1.589 

0.384 

1.540 

0.388 

1.504 

0.391 

1.492 

0.392 

25.0000 

1.621 

0.363 

1.580 

0.360 

1.551 

0.357 

1.541 

0.357 

30.0000 

1.626 

0.351 

1.592 

0.341 

1.568 

0.333 

1.559 

0.336 

35.0000 

1.627 

0.373 

1.584 

0.356 

1.553 

0.344 

1.542 

0.348 

40.0000 

1.622 

0.428 

1.575 

0.409 

1.542 

0.395 

1.530 

0.398 
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Table  A-5.  Index  of  refraction  as  a  function  of  relative  humidity  (0-80  %)  for 
small  rural  aerosols 


Relative  Humiditvf%)  1 

A 

0 

50 

70 

80 

n 

k 

n 

k 

n 

k 

n 

k 

0.3371 

1.530 

0.00590 

1.520 

0.00560 

1.503 

0.00504 

1.449 

0.00331 

0.4000 

1.530 

0.00590 

1.520 

0.00560 

1.502 

0.00504 

1.446 

0.00331 

0.4880 

1.530 

0.00590 

1.520 

0.00560 

1.501 

0.00504 

1.444 

0.00331 

0.5145 

1.530 

0.00590 

1.520 

0.00560 

1.501 

0.00504 

1.444 

0.00331 

0.5500 

1.530 

0.00660 

1.520 

0.00626 

1.501 

0.00563 

1.443 

0.00370 

0.6328 

1.530 

0.00660 

1.520 

0.00626 

1.501 

0.00563 

1.443 

0.00370 

0.6943 

1.530 

0.00730 

1.520 

0.00692 

1.501 

0.00623 

1.443 

0.00409 

0.8600 

1.520 

0.0108 

1.510 

0.0102 

1.492 

0.00922 

1.436 

0.00606 

1.0600 

1.520 

0.0143 

1.510 

0.0136 

1.492 

0.0122 

1.435 

0.00802 

3.0000 

1.342 

0.0190 

1.343 

0.0320 

1.346 

0.0560 

1.355 

0.130 

3.5000 

1.399 

0.00680 

1.399 

0.00693 

1.399 

0.00718 

1.399 

0.00794 

4.0000 

1.397 

0.00710 

1.394 

0.00697 

1.390 

0.00673 

1.377 

0.00600 

4.5000 

1.400 

0.0133 

1.397 

0.0133 

1.390 

0.0133 

1.370 

0.0133 

5.0000 

1.390 

0.0132 

1.387 

0.0132 

1.380 

0.0131 

1.361 

0.0128 

7.9000 

1.185 

0.0575 

1.191 

0.0563 

1.201 

0.0540 

1.233 

0.0471 

8.2000 

1.046 

0.0922 

1.058 

0.0893 

1.081 

0.0839 

1.151 

0.0671 

8.5000 

1.300 

0.178 

1.299 

0.170 

1.297 

0.157 

1.290 

0.116 

9.0000 

2.302 

0.301 

2.249 

0.288 

2.150 

0.263 

1.845 

0.186 

9.5000 

1.884 

0.161 

1.851 

0.155 

1.790 

0.144 

1.602 

0.110 

10.0000 

1.799 

0.112 

1.769 

0.108 

1.714 

0.103 

1.544 

0.0849 

10.5910 

1.718 

0.0850 

1.690 

0.0841 

1.639 

0.0824 

1.481 

0.0773 

11.0000 

1.690 

0.0665 

1.662 

0.0681 

1.611 

0.0709 

1.454 

0.0798 

11.5000 

1.646 

0.0629 

1.619 

0.0670 

1.570 

0.0745 

1.418 

0.0976 

12.5000 

1.587 

0.0641 

1.563 

0.0741 

1.519 

0.0926 

1.383 

0.150 

14.0000 

1.548 

0.0766 

1.531 

0.0917 

1.499 

0.120 

1.400 

0.205 

15.0000 

1.465 

0.170 

1.455 

0.182 

1.436 

0.204 

1.379 

0.272 

18.0000 

1.878 

0.161 

1.855 

0.174 

1.811 

0.199 

1.678 

0.277 

20.0000 

1.988 

0.220 

1.962 

0.229 

1.914 

0.245 

1.765 

0.296 

25.0000 

1.907 

0.268 

1.888 

0.273 

1.852 

0.281 

1.742 

0.307 

30.0000 

1.814 

0.336 

1.800 

0.336 

1.776 

0.335 

1.698 

0.332 

35.0000 

1.914 

0.430 

1.894 

0.425 

1.858 

0.416 

1.746 

0.389 

40.0000 

1.932 

0.530 

1.911 

0.523 

1.872 

0.509 

1.751 

0.466 

Table  A-6.  Index  of  refraction  as  a  function  of  relative  humidity  (90-99  %)  for 
small  rural  aerosols 


■■ 

Relative  Humiditv(’%')  I 

90 

95 

98 

99  1 

Emm 

n 

k 

n 

k 

n 

k 

n 

k 

0.3371 

1.407 

0.00198 

1.393 

0.00153 

1.379 

0.00108 

1.371 

0.000819 

0.4000 

1.403 

0.00198 

1.388 

0.00153 

1.374 

0.00108 

1.366 

0.000819 

0.4880 

1.401 

0.00198 

1.385 

0.00153 

1.371 

0.00108 

1.362 

0.000819 

0.5145 

1.400 

0.00198 

1.385 

0.00153 

1.370 

0.00108 

1.361 

0.000819 

0.5500 

1.399 

0.00222 

1.384 

0.00171 

1.369 

0.00121 

1.360 

0.000916 

0.6328 

1.399 

0.00222 

1.383 

0.00171 

1.368 

0.00121 

1.359 

0.000916 

0.6943 

1.398 

0.00245 

1.382 

0.00189 

1.368 

0.00134 

1.359 

0.00101 

0.8600 

1.393 

0.00363 

1.378 

0.00279 

1.364 

0.00198 

1.356 

0.00150 

1.0600 

1.391 

0.00481 

1.376 

0.00370 

1.362 

0.00263 

1.353 

0.00199 

1.361 

0.187 

1.364 

0.207 

1.366 

0.266 

1.367 

0.237 

1.400 

0.00853 

1.400 

0.00873 

1.400 

0.00892 

1.400 

0.00904 

1.366 

0.00544 

1.363 

0.00525 

1.359 

1.357 

0.00495 

1.355 

0.0134 

1.350 

0.0134 

1.344 

1.341 

0.0134 

1.347 

0.0127 

1.342 

0.0126 

1.337 

1.344 

0.0125 

1.257 

0.0418 

1.266 

0.0400 

1.274 

1.279 

0.0372 

1.205 

0.0543 

1.224 

0.0499 

1.242 

1.253 

0.0430 

1.285 

0.0840 

1.284 

0.0731 

1.282 

0.0625 

1.281 

0.0562 

■aBBiiBI 

1.611 

0.128 

1.531 

0.107 

1.453 

0.0878 

1.406 

0.0761 

9.5000 

1.458 

0.0834 

1.409 

0.0744 

1.361 

0.0657 

1.332 

0.0605 

10.0000 

1.413 

0.0712 

1.368 

0.0665 

1.325 

0.0620 

1.299 

0.0592 

10.5910 

1.360 

0.0733 

1.318 

0.0720 

1.278 

0.0706 

1.254 

0.0699 

11.0000 

1.333 

0.0866 

1.292 

0.0890 

1.252 

0.0912 

1.228 

0.0926 

11.5000 

1.301 

0.115 

1.260 

0.122 

1.221 

0.127 

1.198 

0.131 

12.5000 

1.279 

0.194 

1.243 

0.209 

1.208 

0.223 

1.187 

0.232 

14.0000 

1.324 

0.271 

1.297 

0.294 

1.272 

0.316 

1.257 

0.329 

15.0000 

1.336 

0.324 

1.320 

0.342 

1.306 

0.359 

1.297 

0.370 

18.0000 

1.576 

0.337 

1.541 

0.357 

1.507 

0.377 

1.486 

0.389 

20.0000 

1.651 

0.335 

1.611 

0.348 

1.573 

0.361 

1.550 

0.369 

25.0000 

1.657 

0.326 

1.628 

0.333 

1.600 

0.340 

1.583 

0.344 

30.0000 

1.639 

0.331 

1.619 

0.330 

1.599 

0.329 

1.587 

0.329 

35.0000 

1.660 

0.368 

1.631 

0.360 

1.602 

0.353 

1.585 

0.349 

40.0000 

1.658 

0.434 

1.626 

0.422 

1.595 

0.412 

1.576 

0.405 
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Table  A-7.  Index  of  refraction  as  a  function  of  relative  humidity  (0-80  %)  for 
large  rural  aerosols 


Relative  Humiditvf%)  1 

A 

0 

50 

70 

80  1 

..  (/^m) 

n 

k 

n 

k 

n 

k 

n 

k 

0.3371 

1.53 

0.00590 

1.520 

0.00559 

1.499 

0.00491 

1.435 

0.00286 

0.4000 

1.53 

0.00590 

1.520 

0.00559 

1.498 

0.00491 

1.431 

0.00286 

0.4880 

1.53 

0.00590 

1.520 

0.00559 

1.497 

0.00491 

1.429 

0.00286 

0.5145 

1.53 

0.00590 

1.520 

0.00559 

1.497 

0.00491 

1.429 

0.00286 

0.5500 

1.53 

0.00660 

1.520 

0.00626 

1.497 

0.00549 

1.428 

0.00319 

0.6328 

1.53 

0.00660 

1.520 

0.00626 

1.497 

0.00549 

1.428 

0.00319 

0.6943 

1.53 

0.00730 

1.520 

0.00692 

1.497 

0.00608 

1.427 

0.00353 

0.8600 

1.52 

0.0108 

1.510 

0.0102 

1.488 

0.00899 

1.421 

0.00523 

1.0600 

1.52 

0.0143 

1.510 

0.0136 

1.487 

0.0119 

1.420 

0.00692 

3.0000 

1.342 

0.0190 

1.344 

0.0321 

1.347 

0.0614 

1.357 

0.150 

3.5000 

1.399 

0.00680 

1.399 

0.00693 

1.399 

0.00724 

1.400 

0.00814 

4.0000 

1.379 

0.00710 

1.394 

0.00697 

1.389 

0.00668 

1.373 

0.00581 

4.5000 

1.40 

0.0133 

1.396 

0.0133 

1.389 

0.0133 

1.365 

0.0134 

5.0000 

1.39 

0.0132 

1.387 

0.0132 

1.397 

0.0131 

1.356 

0.0128 

7.9000 

1.185 

0.0575 

1.191 

0.0563 

1.203 

0.0535 

1.241 

0.0453 

8.2000 

1.046 

0.0922 

1.058 

0.0892 

1.086 

0.0826 

1.170 

0.0627 

8.5000 

1.30 

0.178 

1.299 

0.170 

1.296 

0.154 

1.289 

0.105 

9.0000 

2.302 

0.301 

2.248 

0.287 

2.128 

0.257 

1.765 

0.166 

9.5000 

1.884 

0.161 

1.851 

0.155 

1.777 

0.141 

1.553 

0.101 

10.0000 

1.799 

0.112 

1.769 

0.108 

1.702 

0.101 

1.499 

0.0802 

10.5910 

1.718 

0.0850 

1.690 

0.0841 

1.628 

0.0821 

1.440 

0.0759 

11.0000 

1.69 

0.0665 

1.662 

0.0681 

1.600 

0.0716 

1.413 

0.0821 

11.5000 

1.646 

0.0629 

1.619 

0.0670 

1.559 

0.0762 

1.378 

0.104 

12.5000 

1.587 

0.0641 

1.563 

0.0742 

1.509 

0.0968 

1.348 

0.165 

14.0000 

1.548 

0.0766 

1.531 

0.0918 

1.491 

0.126 

1.374 

0.228 

15.0000 

1.465 

0.170 

1.455 

0.182 

1.432 

0.209 

1.364 

0.290 

18.0000 

1.878 

0.161 

1.854 

0.174 

1.802 

0.205 

1.643 

0.298 

20.0000 

1.988 

0.220 

1.962 

0.228 

1.903 

0.249 

1.726 

0.309 

25.0000 

1.907 

0.268 

1.888 

0.273 

1.844 

0.283 

1.713 

0.313 

30.0000 

1.814 

0.336 

1.800 

0.336 

1.770 

0.335 

1.678 

0.332 

35.0000 

1.914 

0.430 

1.894 

0.425 

1.850 

0.414 

1.717 

0.381 

40.0000 

1.932 

0.530 

1.911 

0.522 

1.863 

0.506 

1.719 

0.455 

Table  A-8.  Index  of  refraction  as  a  function  of  relative  humidity  (90-99  %)  for 
large  rural  aerosols 


Relative  Humiditvf%)  1 

A 

90 

95 

98 

99  1 

(Min) 

n 

k 

n 

k 

n 

k 

n 

k 

0.3371 

1.40 

0.00174 

1.386 

0.00132 

1.361 

0.000509 

1.354 

0.000289 

0.4000 

1.395 

0.00174 

1.382 

0.00132 

1.355 

0.000509 

1.348 

0.000289 

0.4880 

1.392 

0.00174 

1.379 

0.00132 

1.352 

0.000509 

1.345 

0.000289 

0.5145 

1.392 

0.00174 

1.378 

0.00132 

1.351 

0.000509 

1.344 

0.000289 

0.5500 

1.391 

0.00194 

1.377 

0.00148 

1.350 

0.000570 

1.343 

0.000323 

0.6328 

1.39 

0.00194 

1.376 

0.00148 

1.349 

0.000570 

1.342 

0.000323 

0.6943 

1.39 

0.00215 

1.376 

0.00164 

1.348 

0.000630 

1.341 

0.000357 

0.8600 

1.385 

0.00318 

1.372 

0.00242 

1.345 

0.000933 

1.338 

0.000529 

1.0600 

1.383 

0.00422 

1.369 

0.00321 

1.343 

0.100 

1.335 

0.000240 

3.0000 

1.362 

0.197 

1.364 

0.215 

1.368 

0.250 

1.370 

0.260 

3.5000 

1.40 

0.00863 

1.400 

0.00882 

1.400 

0.00918 

1.400 

0.00927 

4.0000 

1.364 

0.00534 

1.361 

0.00516 

1.355 

0.00482 

1.353 

0.00472 

4.5000 

1.352 

0.0134 

1.347 

0.0134 

1.338 

0.0134 

1.335 

0.0134 

5.0000 

1.344 

0.0126 

1.340 

0.0126 

1.331 

0.0126 

1.328 

0.0125 

7.9000 

1.262 

0.0409 

1.270 

0.0392 

1.285 

0.0359 

1.289 

0.0351 

8.2000 

1.215 

0.0519 

1.232 

0.0479 

1.265 

0.0400 

1.274 

0.0379 

8.5000 

1.284 

0.0781 

1.283 

0.0682 

1.280 

0.0488 

1.279 

0.0435 

9.0000 

1.568 

0.117 

1.495 

0.0984 

1.352 

0.0624 

1.313 

0.0527 

9.5000 

1.432 

0.0786 

1.387 

0.0704 

1.298 

0.0544 

1.274 

0.0500 

10.0000 

1.389 

0.0687 

1.348 

0.0644 

1.268 

0.0561 

1.246 

0.0538 

10.5910 

1.338 

0.0726 

1.300 

0.0714 

1.226 

0.0689 

1.205 

0.0683 

11.0000 

1.311 

0.0879 

1.273 

0.0900 

1.199 

0.0942 

1.179 

0.0953 

11.5000 

1.279 

0.119 

1.243 

0.124 

1.171 

0.135 

1.151 

0.138 

12.5000 

1.26 

0.202 

1.227 

0.215 

1.163 

0.242 

1.146 

0.249 

14.0000 

1.31 

0.284 

1.286 

0.304 

1.239 

0.345 

1.227 

0.356 

15.0000 

1.327 

0.334 

1.314 

0.350 

1.287 

0.382 

1.280 

0.391 

18.0000 

1.557 

0.348 

1.525 

0.366 

1.462 

0.403 

1.445 

0.413 

20.0000 

1.63 

0.342 

1.594 

0.354 

1.524 

0.378 

1.505 

0.385 

25.0000 

1.642 

0.330 

1.615 

0.336 

1.563 

0.348 

1.549 

0.352 

30.0000 

1.628 

0.330 

1.610 

0.330 

1.574 

0.329 

1.564 

0.328 

35.0000 

1.645 

0.364 

1.618 

0.357 

1.565 

0.344 

1.551 

0.341 

40.0000 

1.641 

0.428 

1.612 

0.418 

1.555 

0.398 

1.539 

0.392 
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Table  A-9.  Index  of  refraction  as  a  function  of  relative  humidity  (0-80  %)  for 
oceanic  aerosols 


0.3371 

0.4000 

0.4880 

0.5145 

0.5500 

0.6328 

0.6943 

0.8600 


4.5000 

5.0000 

7.9000 

8.2000 


0 

9.5000 
10.0000 
10.5910 
11.0000 

11.5000 

12.5000 


14.0000 

15.0000 

18.0000 

20.0000 

25.0000 

30.0000 

35.0000 

40.0000 


0 

n 

k 

1.510 

4.00  X  10 

1.500 

3.00  X  10 

1.500 

2.00  X  10 

1.500 

1.00  X  10 

1.500 

1.00  X  10 

1.490 

2.00  X  10 

1.490 

1.00  X  10 

1.480 

3.00  X  10 

Relative  Humiditvf% 


1.470 

1.610 

1.480 

1.480 

1.490 

1.470 

1.400 

1.420 


1.480 

1.650 

1.580 

1.548 

1.500 

1.480 

1.480 

1.420 


1.410 

1.450 

1.780 

1.760 

1.760 

1.770 

1.760 

1.740 


2.00  X  10 
1.00  X  10 
1.60  X  10 
1.40  X  10 


2.60  X  10-2 

2.80  X  10-2 

1.80  X  10-2 
1.50  X  10-2 
1.40  X  10-2 
1.40  X  10-2 
1.40  X  10-2 

1.60  X  10-2 


2.30  X  10-2 
3.50  X  10-2 

1.30  X  10-^ 
1.52  X  10-^ 
2.05  X  10-^ 
3.00  X  10-1 
5.00  X  10-1 

1.00 


..444 

.567 

.465 

.457 

.461 

.444 

.381 

1.396 


1.443 

1.580 

1.519 

1.482 

1.442 

1.421 

1.416 

1.366 


1.374 

1.417 

1.715 

1.709 

1.718 
1.730 

1.719 
1.700 


3.29  X  10-^  1.425  1.97 
2.49  X  10-®  1.417  1.54 
1.65  X  10-®  1.415  1.01 
8.40  X  10-®  1.414  5.43 
8.54x10-^  1.413  5.83 
1.408  1.72 


1.425  1.97  X  10-^ 
1.417  1.54  X  10-® 
1.415  1.01  X  10-® 
1.414  5.43  X  10-^ 
1.413  5.83  X  10-^ 
1.408  1.72  X  10-® 
1.408  6.40  X  10-® 
1.402  1.62  X  10"® 


395  9.85 
486  1.46 


1.90  X  10-® 
8.74  X  10-® 
2.52  X  10-® 


1.64  X  lO-'^ 
5.76  X  10-2 
3.02  X  10-^ 
1.98  X  10-3 
3.58  X  10-3 
4.30  X  10-3 
1.68  X  10-2 
2.27  X  10-2 


2.79  X  10-2  1.375  3.15 
3.02  x  10-2  1.449  3.42 
2.28  x  10-2  1.405  3.17 
2.15  x  10-2  1.373  3.35 
2.37  x  10-2  1.334  4.17 


2.90  X  10-2 
3.72  X  10-2 
6.01  X  10-2 


8.60  X  10-2 
1.02  X  10-1 
1.84  X  10-1 
1.96  X  10-1 
2.32  X  10-1 
3.05  X  10-1 
4.70  X  10-1 

8.88  X  10-1 


1.311  5.69 
1.297  8.03 
1.266  1.42 


1.306  2.03 
1.357  2.25 
1.595  2.83 
1.615  2.77 

1.641  2.83 
1.657  3.15 

1.642  4.15 
1.626  6.81 


X  10  2 

X  10-2 

X  10-2 
X  10-2 
X  10-2 
X  10-2 

X  10-2 

X  10-^ 


X  10-^ 
X  10-^ 
X  10-^ 
X  10-3 

X  10-3 
X  10-3 


n 


1.366 

1.359 

1.356 

1.355 

1.354 

1.352 

1.351 

1.348 


5.83  X  10-® 
5.44  X  10-® 
3.39  X  10-^ 
2.30  X  10-® 
2.98  X  10-® 
1.53  X  10-® 
3.94  X  10-® 
6.69  X  10 


L0-® 

1.344 

2.91  X  10-® 

10-3 

1.401 

2.39  X  10-3 

L0-® 

1.410 

8.41  X  10-® 

CO 

1 

O 

1.367 

4.19  X  10-® 

L0-® 

1.352 

1.19  X  10-2 

[0-® 

1.343 

1.11  X  10-2 

O 

1 

to 

1.307 

3.12  X  10-2 

.0-2 

1.303 

3.32  X  10-2 

1.304 

1.311 

1.286 

1.259 

1.220 

1.195 

1.171 

1.161 


1.235 

1.293 

1.468 

1.516 

1.560 
1.579 

1.561 
1.547 


3.53  X  10-2 
3.84  X  10-2 
4.10  X  10-2 

4.62  X  10-2 
6.06  X  10-2 

8.63  X  10-2 
1.26  X  10-3 
2.28  X  10-3 


3.26  X  10-3 
3.55  X  10-3 
3.88  X  10-3 

3.62  X  10-3 
3.37  X  10-3 
3.24  X  10-3 
3.57  X  10-3 

4.63  X  10-3 


Table  A- 10.  Index  of  refraction  as  a  function  of  relative  humidity  (90-99  %)  for 
oceanic  aerosols 


Relative  Humiditv(%)  I 

HI 

90 

95 

98 

99  1 

n 

k 

n 

k 

n 

k 

n 

k 

0.3371 

1.357 

3.76  X  10-® 

1.352 

2.49  X  10-5 

1.348 

1.58  X  10-5 

1.347 

1.22  X  10-5 

0.4000 

1.351 

3.96  X  10-® 

1.346 

3.04  X  10-9 

1.342 

2.39  X  10-9 

1.341 

2.13  X  10-9 

0.4880 

1.347 

2.39  X  10"® 

1.342 

1.77  X  10-9 

1.338 

1.33  X  10-9 

1.337 

1.15  X  10-9 

0.5145 

1.346 

1.83  X  10-® 

1.341 

1.54  X  10-9 

1.337 

1.34  X  10-9 

1.336 

1.26  X  10-9 

0.5500 

1.345 

2.30  X  10-9 

1.336 

2.11  X  10-9 

1.335 

2.04  X  10-9 

0.6328 

1.344 

1.339 

1.49  X  10-5 

1.335 

1.47  X  10-5 

1.334 

1.47  X  10-5 

0.6943 

1.343 

1.338 

3.34  X  10-5 

1.334 

3.18  X  10-5 

1.333 

3.12  X  10-5 

0.8600 

1.335 

4.41  X  10-2 

1.332 

3.79  X  10-^ 

1.330 

3.55  X  10-^ 

1.0600 

1.337 

1.88  X  10-5 

1.332 

1.24  X  10-5 

1.329 

7.85  X  10-5 

1.327 

6.08  X  10-5 

3.0000 

1.389 

2.52  X  10-1 

1.381 

2.61  X  10-^ 

1.375 

2.67  X  10-^ 

1.373 

2.69  X  10-5 

3.5000 

1.406 

8.82  X  10-5 

1.403 

9.07  X  10-5 

1.401 

9.25  X  10-5 

1.401 

9.32  X  10-5 

4.0000 

1.361 

4.36  X  10-5 

1.356 

4.47  X  10-5 

1.353 

4.54  X  10-5 

1.352 

4.57  X  10-5 

4.5000 

1.344 

1.25  X  10-2 

1.339 

1.29  X  10-2 

1.335 

1.32  X  10-2 

1.334 

1.33  X  10-2 

5.0000 

1.336 

1.17  X  10-2 

1.331 

1.20  X  10-2 

1.328 

1.22  X  10-2 

1.326 

1.23  X  10-2 

7.9000 

1.302 

3.23  X  10-2 

1.298 

3.30  X  10-2 

1.296 

3.35  X  10-2 

1.295 

3.37  X  10-2 

8.2000 

1.296 

3.40  X  10-2 

1.292 

3.45  X  10-2 

1.289 

3.48  X  10-2 

1.287 

3.50  X  10-2 

8.5000 

1.293 

3.59  X  10-2 

1.286 

3.62  X  10-2 

1.282 

3.65  X  10-2 

1.280 

3.65  X  10-2 

9.0000 

1.291 

3.90  X  10-2 

1.278 

3.94  X  10-2 

1.269 

3.97  X  10-2 

1.266 

3.98  X  10-2 

9.5000 

1.268 

4.24  X  10-2 

1.257 

4.32  X  10-2 

1.249 

4.39  X  10-2 

1.246 

4.41  X  10-2 

10.0000 

1.242 

4.81  X  10-2 

1.231 

4.93  X  10-2 

1.224 

5.01  X  10-2 

1.221 

5.05  X  10-2 

10.5910 

1.203 

6.34  X  10-2 

1.192 

6.52  X  10-2 

1.185 

6.64  X  10-2 

1.182 

6.69  X  10-2 

11.0000 

1.177 

9.06  X  10-2 

1.167 

9.33  X  10-2 

1.159 

9.52  X  10-2 

1.156 

9.60  X  10-2 

11.5000 

1.152 

1.32  X  10-1 

1.141 

1.37  X  10-1 

1.133 

1.40  X  10"^ 

1.129 

1.41  X  10-5 

12.5000 

1.145 

2.41  X  10-1 

1.135 

2.49  X  10-1 

1.129 

2.54  X  10-5 

1.126 

2.57  X  10-5 

14.0000 

1.225 

3.44  X  10-1 

1.218 

3.55  X  10-1 

1.214 

3.63  X  10-5 

1.212 

3.67  X  10-5 

15.0000 

1.283 

3.75  X  10-1 

1.278 

3.87  X  10-1 

1.273 

3.95  X  10-5 

1.272 

3.98  X  10-5 

18.0000 

1.450 

4.04  X  10-1 

1.438 

4.14  X  10-1 

1.430 

4.20  X  10-5 

1.426 

4.23  X  10-5 

20.0000 

1.501 

3.75  X  10-1 

1.492 

3.83  X  10-1 

1.485 

3.88  X  10-5 

1.483 

3.91  X  10-5 

25.0000 

1.548 

3.45  X  10-1 

1.541 

3.50  X  10-1 

1.535 

3.53  X  10-5 

1.533 

3.55  X  10-5 

30.0000 

1.567 

3.26  X  10-1 

1.560 

3.27  X  10-1 

1.555 

3.27  X  10-5 

1.553 

3.28  X  10-5 

35.0000 

1.549 

3.48  X  10-1 

1.542 

3.43  X  10-1 

1.536 

3.39  X  10-5 

1.534 

3.38  X  10-5 

40.0000 

1.535 

4.31  X  10-1 

1.528 

4.11  X  10-1 

1.523 

3.97  X  10-5 

1.521 

3.91  X  10-5 

57 


This  page  left  intentionally  blank. 


Appendix  B 
GRAPHS 


The  following  pages  contain  three-dimensional  plots  of  the  Phase  FuNction 
DATabase  (PFNDAT)  phase  functions  (on  the  ordinate  axis)  versus  angle  and 
wavelength.  PFNDAT  cases  1  to  30  have  been  evaluated  at  32  different 
wavelengths  from  0.35  to  40  /im  and  the  remainder  of  the  cases  at  16  different 
wavelengths  ranging  from  0.55  to  12  ^m.  The  wavelength  scale  used  in  the  plots 
(the  y  dimension)  is  log-based  as  is  the  vertical  scale  (the  z  dimension)  that 
represents  the  phase  function  value.  In  a  majority  of  the  plots  the  variation  of 
the  zero  peak  phase  function  value  is  linear  with  wavelength.  This  is  due  to 
the  direct  relationship  between  the  amount  of  forward  scatter  and  the  particle 
size  parameter.  The  plots  are  created  using  exactly  the  wavelength  and  angle 
information  contained  in  the  database.  Thus,  where  the  plot  lines  are  sparse,  so 
are  the  computed  data. 
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o 


Figure  B-1.  Maritime  aerosol,  0  percent  relative  humidity. 


Wavelength  (microns) 


Maritime  Aerosol,  50%  Rh  EOSAEL  02 


o 


Figure  B-2.  Maritime  aerosol,  50  percent  relative  humidity. 


Wavelength  (microns) 


o 


Figure  B-3.  Maritime  aerosol,  70  percent  relative  humidity. 


Wavelength  (microns) 


Maritime  Aerosol,  80%  Rh  EOSAEL  04 


o 


Figure  B-4.  Maritime  aerosol,  80  percent  relative  humidity. 


Wavelength  (microns) 


o 


Figure  B-5.  Maritime  aerosol,  90  percent  relative  humidity. 
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Wave 1 eng  th  (microns) 


Maritime  Aerosol,  95%  Rh  EOSAEL  06 


o 


Figure  B-6.  Maritime  aerosol,  95  percent  relative  humidity. 
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Figure  B-7.  Maritime  aerosol,  98  percent  relative  humidity. 
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Wavelength  (microns) 


Maritime  Aerosol,  99%  Rh  EOSAEL  08 


o 

rj* 


Figure  B-8.  Maritime  aerosol,  99  percent  relative  humidity. 


microns 


o 


Figure  B-9.  Urban  aerosol,  0  percent  relative  humidity. 
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Wavelength  (microns) 


o 


Figure  B-10.  Urban  aerosol,  50  percent  relative  humidity. 


Wavelength  (microns) 


o 
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Wavelength  (microns) 


o 


Figure  B-12.  Urban  aerosol,  80  percent  relative  humidity. 


Wavelength  (microns) 


o 
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Wavelength  (microns) 


Urban  Aerosol,  95%  Rh  EOSAEL  14 


o 


Figure  B-14.  Urban  aerosol,  95  percent  relative  humidity. 
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Wavelength  (microns) 


o 


Figure  B-15.  Urban  aerosol,  98  percent  relative  humidity. 
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Wavelength  (microns) 


Urban  Aerosol,  99%  Rh  EOSAEL  16 


o 


Figure  B-16.  Urban  aerosol,  99  percent  relative  humidity. 


Wavelength  (microns) 


o 


Figure  B-17.  Rural  aerosol,  0  percent  relative  humidity. 
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Wavelength  (microns) 


o 


Figure  B-18.  Rural  aerosol,  50  percent  relative  humidity, 


Wavelength  (microns) 


o 


Figure  B-19.  Rural  aerosol,  70  percent  relative  humidity. 
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Figure  B-20.  Rural  aerosol,  80  percent  relative  humidity. 
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Wavelength  (microns) 


o 


Figure  B-21.  Rural  aerosol,  90  percent  relative  humidity. 
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Wavelength  (microns) 


o 
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Wavelength  (microns) 


o 
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Wavelength  (microns) 


o 


Figure  B-24.  Rural  aerosol,  99  percent  relative  humidity. 


Wavelength  (microns) 
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Wavelength  (microns) 
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Figure  B-26.  Fog  (moderate  radiation). 
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Wavelength  (microns) 


Drizzle,  Inmi/hr  EOSAEL  27 


O 


Figure  B-27.  Rain  (drizzle) 


Wavelength  (microns) 


widespread  Rain,  5iran/hr  EOSAEL  28 


o 


Figure  B-28.  Rain  (widespread) 


o 


Figure  B-29.  Rain  (thunderstorm). 
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Wavelength  (microns) 


o 


Figure  B-30.  Snow 


Wavelength  (microns) 


o 


Figure  B-31.  Dust  (light  aerosol  loading) 


Wavelength  (microns) 


Heavy  Loading  Dust  EOSAEL  51 


o 


Figure  B-32.  Dust  (heavy  aerosol  loading). 


o 


Figure  B-34.  Smoke  (white  phosphorous),  17  percent  relative  humidity. 
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Wavelength  (microns) 


o 


Figure  B-35.  Smoke  (white  phosphorous),  50  percent  relative  humidity. 


Wavelength  (microns) 


o 


in 

in 

W 

< 


CO 

O 

W 


I 


Figure  B-36,  Smoke  (white  phosphorous),  90  percent  relative  humidity. 


Wavelength  (microns) 


Figure  B-37.  Smoke  (fog  oil),  50  percent  relative  humidity. 


Wavelength  (microns) 


HC  Smoke,  85%  Rh  EOSAEL  57 


o 


Figure  B-38.  Smoke  (hexachloroethane),  85  percent  relative  humidity. 


Wavelength  (microns) 
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Appendix  C 
EOSAEL  MODULES 


AGAUS 
BITS 
AS  CAT 
CL  I  MAT 
CLTRAN 
COMBIC 

COPTER 

FASCAT 

FCLOUD 

OVRCST 

FITTE 

GRNADE 

GSCAT 

ILUMA 

IMTURB 

KWIK 

LASS 

LOWTRIM 

LZTRAN 

MPLUME 

NBSCAT 

MSCAT 

NMMW 

NOVAE 

OVRCST 

PFNDAT 

RADAR 

REFRAC 

TARGAC 

UVTRAN 

XSCALE 


Mie  Scattering  Code 
Broad-band  Integrated  Transmittances 
Approximate  Multiple  Scattering 
Climatology 

Transmission  through  Clouds 

Obscuration  Model  for  Multiple  Battlefield-Induced 
Contaminants 

Obscuration  due  to  Helicopter-Lofted  Snow  and  Dust 
Fast  Algorithm  for  Atmospheric  Scattering  Calculations 
Transmission  through  Cloud  of  Ellipsoidal  Geometry 
Path  Radiance/Contrast  Beneath  Overcast  Conditions 
Fire- Induced  Transmission  and  Turbulence  Effects 
Smoke  Munitions  Self-Screening  Applications 
Multiple  Scattering  using  Gaussian  Geometry 
Natural  Illumination  under  Realistic  Weather  Conditions 
Imaging  Through  Optical  Turbulence 
Transmission  Threshold  Smoke  Munitions  Expenditures 
Model 

Large  Area  Screening  Systems  Application 
Atmospheric  Transmittance  and  Radiance  for 
Broadband  Applications 

Laser  Transmittance-Gaseous  Absorption  Algorithm 

Missile  Smoke  Plume  Obscuration 

Narrow  Beam  Multiple  Scattering 

Aerosol  Multiple  Scattering,  Monte  Carlo 

Near  Millimeter  Wave,  Gaseous  Absorption 

Nonlinear  Aerosol  Vaporization  and  Breakdown  Effects, 

High  Energy  Lasers 

Contrast  Transmission 

Aerosol  Phase  Function  Data  Base 

Millimeter  Wave  System  Performance 

Optical  Path  Bending  Code  for  Near  Earth  Paths 

Target  Acquisition 

Ultraviolet  Transmission  and  Lidar  Simulation 
Natural  Aerosol  Extinction 
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Distribution 


Copies 


AEROSPACE 

ATTN  M  DICKERSON  1 

PO  BOX  92957 

LOS  ANGELES  CA  90009  2957 

AEROSPACE 

M4/954 

ATTNC  YINGER  1 

PO  BOX  92957 

LOS  ANGELES  CA  90009 

AEROSPACE  CORP 
M4/954 

ATTN  J  LANGER  1 

PO  BOX  92957 

LOS  ANGELES  CA  90009 

AEROSPACE  CORP 
MSM4944 

ATTN  JTOTH  1 

PO  BOX  92957 

LOS  ANGELES  CA  90009  2957 

AFMC  DOW  1 

WRIGHT  PATTERSON  AFB  OH  0334  5000 

AFSPC/DRFN 

ATTNCAPTRKOON  1 

150  VANDENBERG  STRET  STE  1 105 
WRIGHT  PATTERSON  AFB  CO  80914  45900 

AIR  WEATHER  SERVICE  1 

TECH  LIBRARY  FL4414  3 
SCOTT  AFB  IL  62225  5458 
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ALLEN’S  TIME 
ATTN  D  ALLAN 
POBOX  66 

FOUNTAIN  GREEN  UT  84632 

AMCOMP/746TS 
ATTNACHASKO 
POBOX  3 10 
fflGH  ROLLS  NM  88325 

APPLIED  RSRCH  LAB  UNIV  OF  TX 
ATTOB  RENFRO 
POBOX  8029 
AUSTIN  TX  78713  8029 

APPLIED  RSRCH  LAB  UNIV  OF  TX 
ATTN  J  SAUNDERS 
POBOX  8029 
AUSTIN  TX  78713  8029 

APPLIED  RSRCH  LAB  UNIV  OF  TX 
ATTN  R  MACH 
POBOX  8029 
AUSTIN  TX  78713  8029 

ARINC 

ATTN  L  CONOVER 
1 925  AEROTECH  DR  STE  2 1 2 
COLORADO  SPRINGS  CO  809 1 6 

ARINC 

ATTN  P  MENDOZA 
4055  HANCOCK  STRET 
SAN  DIEGO  CA  921 10 

ARL  CHEMICAL  BIOLOGY 
NUC  EFFECTS  DIV 
AMSRL  SL  CO 
APGMD  21010  5423 
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ARMY  ARDEC 
SMCAR IMI I  BLDG  59 
DOVER  NJ  07806  5000 

ARMY  COMMUNICATIONS 
ELECTR  CTR  FOR  EW  RSTA 
AMSEL  EW  D 

FT  MONMOUTH  NJ  07703  5303 

ARMY  COMMUNICATIONS 
ELECTR  CTR  FOR  EW  RSTA 
AMSEL  EW  MD 
FT  MONMOUTH  NJ  07703  5303 

ARMY  CORPS  OF  ENGRS 
ENGR  TOPOGRAPHICS  LAB 
ETL  GS  LB 

FT  BELVOIRVA  22060 

ARMY  DUGWAY  PROVING  GRD 
STEDP3 

DUGWAY  UT  84022  5000 

ARMY  DUGWAY  PROVING  GRD 
STEDP  MT  DA  L  3 
DUGWAY  UT  84022  5000 

ARMY  DUGWAY  PROVING  GRD 
STEDP  MTM 
ATTN  MR  BOWERS 
DUGWAY  UT  84022  5000 

ARMY  FIELD  ARTLLRY  SCHOOL 
ATSF  TSM  TA 
FT  SILL  OK  73503  5600 

ARMY  FOREGN  SCI  TECH  CTR 
CM 

220  7TH  STREET  NE 
CHARLOTTESVILLE  VA  22901  5396 


ARMY  INFANTRY 
ATSH  CD  CS  OR 
ATTN  DR  E  DUTOIT 
FT  BENNING  GA  30905  5090 

ARMY  MATERIEL  SYST 
ANALYSIS  ACTIVITY 
AMXSY 

APGMD  21005  5071 

ARMY  MATERIEL  SYST 
ANALYSIS  ACTIVITY 
AMXSY  AT 
ATTN  MR  CAMPBELL 
APG  MD  21005  5071 

ARMY  MATERIEL  SYST 
ANALYSIS  ACTIVITY 
AMXSY  CR 
ATTNMRMARCHET 
APGMD 21005  5071 

ARMY  MATERIEL  SYST 
ANALYSIS  ACTIVITY 
AMXSY  CS 
ATTN  MR  BRADLEY 
APGMD 21005  5071 

ARMY  MISSILE  CMND 
AMSMI 

REDSTONE  ARSENAL  AL  35898  5243 

ARMY  MISSILE  CMND 
AMSMI 

REDSTONE  ARSENAL  AL  35898  5253 

ARMY  MISSILE  CMND 
AMSMI  RD  AC  AD 
ATTN  DR  PETERSON 
REDSTONE  ARSENAL  AL  35898  5242 
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ARMY  MISSILE  CMND 
AMSMI  RD  AS  SS 
ATTN  MR  H  F  ANDERSON 
REDSTONE  ARSENAL  AL  35898  5253 

ARMY  MISSILE  CMND 
AMSMI  RD  AS  SS 
ATTN  MR  B  WILLIAMS 
REDSTONE  ARSENAL  AL  35898  5253 

ARMY  MISSILE  CMND 
AMSMI  RD  DE  SE 
ATTN  MR  GORDON  LILL  JR 
REDSTONE  ARSENAL  AL  35898  5245 

ARMY  MISSILE  CMND 
REDSTONE  SCI  INFO  CTR 
AMSMI  RDCSR  DOC 
REDSTONE  ARSENAL  AL  35898  5241 

ARMY  NUCLEAR  CML  AGCY 
MONA  ZB  BLDG  2073 
SPRINGFIELD  VA  22150  3198 

ARMYOEC 

CSTE  EFS 

PARK  CENTER  IV 

4501  FORD  AVE 

ALEXANDRIA  VA  22302  1458 

ARMY  RESEARCH  LABORATORY 
AMSRL 

2800  POWDER  MILL  ROAD 
ADELPHI  MD  20783  1145 

ARMY  RESEARCH  LABORATORY 
AMSRL  BE  E 

BATTLEFIELD  ENVIR  DIR 
WSMRNM  88002  5501 


ARMY  RESEARCH  LABORATORY 
AMSRL  BE  S 

BATTLEFIELD  ENVIR  DIR 
WSMRNM  88002  5501 

ARMY  RESEARCH  LABORATORY 
AMSRL  BE  W 
BATTLEFIELD  ENVIR  DIR 
WSMRNM  88002  5501 

ARMY  RESEARCH  LABORATORY 
AMSRL  D 

2800  POWDER  MILL  ROAD 
ADELPHI  MD  20783  1 145 

ARMY  RESEARCH  LABORATORY 
AMSRL  OP  Cl  SD  TL 
2800  POWDER  MILL  ROAD 
ADELPHI  MD  20783  1145 

ARMY  RESEARCH  LABORATORY 
AMSRL  OP  SD  TP 
TECHNICAL  PUBLISHING 
2800  POWDER  MILL  ROAD 
ADELPHI  MD  20783  1145 

ARMY  RESEARCH  LABORATORY 
AMSRL  SS  SH 
ATTN  DR  SZTANKAY 
2800  POWDER  MILL  ROAD 
ADELPHI  MD  20783  1 145 

ARMY  RESEARCH  LABORATORY 
AMSRL  IS  ES 
ATTN  ALAN  WETMORE 
2800  POWDER  MILL  ROAD 
ADELPHI  MD  20783  1 145 
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1 


ARMY  RESEARCH  LABORATORY 
ELECTROMAG  GROUP 

CAMPUS  MAIL  CODE  F  0250  UNIVERSITY  OF  TX 
ATTN  A  TUCKER 
AUSTIN  TX  78712 

ARMY  RESEARCH  OFC  1 

ATTN  AMXRO  GS 

ATTN  DR  BACH 

PO  BOX  12211 

RTP  NC  27009 

ARMY  SATELLITE  COMM  AGCY  1 

DRCPM  SC  3 

FT  MONMOUTH  NJ  07703  5303 

ARMY  STRAT  DEFNS  CMND  1 

CSSD  SL  L 
ATTN  DR  LILLY 
POBOX  1500 

HUNTSVILLE  AL  35807  3801 

ARMY  TOPO  ENGR  CTR  1 

CETEC  ZC  1 

FT  BELVOIR  VA  22060  5546 

ARMY  TRADOC  ANALYSIS  CTR  1 

ATRC  WSS  R 
WSMRNM  88002  5502 

ARPA 

ATTN  L  STOTTS  1 

JPENNELLA 

B  KASPAR 

3701  N  FAIRFAX  DR 

ARLINGTON  VA  22203  1 7 1 4 
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ASC  OLATJH 
JDAM  PIP 

ATTN  LTV  JOLLEY 
102  W  DAVE 
EGLINAFBFL  32542 

ASHTECHINC 
ATTN  S  GOUREVITCH 
1177KIFERRD 
SUNNYVALE  CA  94086 

ASST  DPTY  CG  FOR  RDE  HDQTRS 
US  ARMY  MATL  CMND 
AMCRD 

ATTN  COL  SMANESS 
5001  EISENHOWER  AVE 
ALEXANDRIA  VA  22333  0001 

ATA  ASSOC 
ATTN  M  HARKINS 
6800BACKLICKRD 
SPRINGFIELD  VA  22 1 50 

ATMOSPHERIC  SCI  DIV 
GEOPHYSICS  DIRCTRT 
PHILLIPS  LABORATORY 
HANSCOM  AFB  MA  01731  5000 

BD  SYSTEMS 
ATTN  J  BUTTS 
385  VAN  NESS  AVE  #200 
TORRANCE  CA  90501 

CECOM 
PM  GPS 

ATTN  COLS  YOUNG 
FT  MONMOUTH  NJ  07703 
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1 


CMD  (420000D(C0245)) 

ATTN  DR  A  SHLANTA 
NAVAIR  WARCENWPNDIV 
1  ADMIN  CIR 

CFIINA  LAKE  CA  93555  6001 


DAMI POI  1 

WASH  DC  20310  1067 

DEAN  RMD  1 

ATTN  DR  GOMEZ 
WASH  DC  20314 

DEPT  OF  COMMERCE  CTR  1 

325  S  BROADWAY 
BOULDER  CO  80303 

DEPT  OF  COMMERCE  CTR  1 


MOUNTAIN  ADMINISTRATION 
SPPRT  CTR  LIBRARY  R  5 1 
325  S  BROADWAY 
BOULDER  CO  80303 

DEPT  OF  MATHEMATICAL  SCIENCES  1 

ATTN  MDN  A  (MAJ  DON  ENGEN) 

THAYER  HALL 

WEST  POINT  NY  10996  1786 

DEPT  OF  THE  AIR  FORCE  1 

OL  A  2D  WEATHER  SQUAD  MAC 
HOLLOMAN  AFB  NM  88330  5000 

DMA 

ATTNRKLOTZ 
4211  BRIARS  RD 
ONLEYMD  20832  1814 


1 


DMA 

LAI 

ATTNBTALLMAN 
3200  2ND  STRET 
ST  LOUIS  MO  63116 

DMA 

L41 

ATTN  B  HAGAN 
3200  S  2ND  STRET 
ST  LOUIS  MO  63118 

DMA 

L41 

ATTN  D  MORGAN,  F  MUELLER 
3200  S  2ND  STRET 
STL0UISM063118 

DOTAFSPC/DRFN 
ATTNHSKALSKI 
1 50  VANDENBERG  STET 
PETERSON  AFB  CO  809 1 4 

DPTY  CG  FOR  RDE  HDQTRS 
US  ARMY  MATL  CMND 
AMCRD 

ATTN  BG  BEAUCHAMP 
5001  EISENHOWER  AVE 
ALEXANDRIA  VA  22333  0001 

DPTY  ASSIST  SCY  FOR  RSRCH  &  TECHL 

SARDTT 

ATTN  DCHAIT 

THE  PENTAGON 

WASHINGTON  DC  203 1 0  01 03 


no 


DPTY  ASSIST  SCY  FOR  RSRCH  &  TECHL 
SARDTT 

ATTN  F  MILTON  RM 
3E479 

THE  PENTAGON 
WASHINGTON  DC  203 1 0  0 1 03 

DPTY  ASSIST  SCY  FOR  RSRCH  &  TECHL 

SARDTT 

ATTN  TRILLION 

THE  PENTAGON 

WASHINGTON  DC  203 1 0  0 1 03 

DPTY  ASSIST  SCY  FOR  RSRCH  &  TECHL 
SARDTT 

ATTNKKOMINOS 
THE  PENTAGON 
WASHINGTON  DC  203 1 0  01 03 

DPTY  ASSIST  SCY  FOR  RSRCH  &  TECHL 
SARDTT 

ATTNBREISMAN 
THE  PENTAGON 
WASHINGTON  DC  203 1 0  0 1 03 

DR  HANS  J  LIEBE 
NTIA  ITS  S  3 
325  S  BROADWAY 
BOULDER  CO  80303 

DR  JERRY  DAVIS 
NCSU 

PO  BOX  8208 
RALEIGH  NC  27650  8208 


DTIC 

8725  JOHN  J.  KINGMAN  RD 
STE  0944 

FT  BELVOIR  VA  22060  62 1 8 

DUSD  SPACE 
1E765 

ATTNJGMCNEFF 

3900  DEFENSE  PENTAGON 

WASHINGTON  DC  20301  3900 

ELECTRONIC  SYSTEMS  DIV  DIR 
CECOMRDEC 
ATTNJNIEMELA 
FT  MONMOUTH  NJ  07703 

G  GIBBONS 

859  WILLIAMETTE  STRT 
EUGENE  OR  9740 16806 

GPS  JOINT  PROG  OFC  DIR 
ATTN  COL  J  CLAY 
2435  VELA  WAY  STE  1613 
LOS  ANGELES  AFB  CA  90245  5500 

HEWLETT  PACKARD  CO 
ATTNJKUSTERS 
5301  STEVENS  CREED  BLVD 
SANTA  CLARA  CA  95052 

HOLLOMAN  AFB 
ATTNKWERNIE 
1 644  V  ANDERGRIFT  RD 
HOLLOMAN  AFB  NM  88330  7850 

HQ  AWS  DOO  1 
SCOTT  AFB  IL  62225  5008 
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HUGHES  AIRCRAFT 
ATTO  R  MALLA 
800  APOLLO  AVE 
PO  BOX  902 
EL  SEGUNDO  CA  90245 

HUGHES  AIRCRAFT 
ATTN  S  PECK 
800  APOLLO  AVE 
PO  BOX  902 
EL  SEGUNDO  CA  90245 

HUGHES  SPACE  &  COMM 

MS/SC/SIO/S364 

ATTN  C  SHEKLELLS 

PO  BOX  92919  AIRPORT  STATION 

LOS  ANGELES  CA  90009 

INTERMETRICS  INC 

ATTN  J  MCGOWAN 

615  HOPE  RD  BLDG  4  2ND  FLOOR 

EATONTOWN  NJ  07724 

ITT  AEROSPACE 
ATTN  P  BRODIE 
100  KINGSLANDRD 
CLIFTON  NJ  07014 

ITT  AEROSPACE 
MS  8538 
ATTN  L  DOYLE 
100  KINGSLANDRD 
CLIFTON  NJ  07014 

ITT  AEROSPACE 
MS  2511 

ATTN  R  TELLER 
100  KINGSLANDRD 
CLIFTON  NM  07014 


ITT  AEROSPACE 
MS  8528 

ATTN  H  RAWICZ 
100  KINGSLAND  RD 
CLIFTON  NJ  07014 

KERNCO 
ATTN  R  KERN 
28  HARBOR  STRET 
DANVERS  MA  01923 

LOCKHEED  MARTIN 
ATTN  B  MARQUIS 
1250  ACADEMY  PARK  LOOP  101 
COLORADO  SPRINGS  CO  80912 

LOCKHEED  MARTIN 
ATTN  J  TAYLOR 

1250  ACADEMY  PARK  LOOP  STE  101 
COLORADO  SPRINGS  CO  80910 

LORAL 

ATTN  B  MATHON 

700  N  FREDERICK  PIKE 

GAITHERSBURG  MD  20879 

LORAL 
B  PISOR 

2915  BASELINE  RD  530 
BOULDER  CO  80303 

LORAL  FEDERAL  SYSTEMS 
ATTN  R  ASTALOS 
9970  FEDERAL  DR 
COLORADO  SPRINGS  CO  80921 

LORAL  FEDERAL  SYSTEMS 
ATTN  M  BAKER 
9970  FEDERAL  DR 
COLORADO  SPRINGS  CO  80921 


LORAL  FEDERAL  SYSTEMS 
ATTN  J  KANE 
9970  FEDERAL  DR 
COLORADO  SPRINGS  CO  80921 

METEOROLOGIST  IN  CHARGE 
KWAJALEIN  MISSILE  RANGE 
PO  BOX  67 

APO  SAN  FRANCISCO  CA  96555 

MIL  ASST  FOR  ENV  SCI  OFC 
OF  THE  UNDERSEC  OF  DEFNS 
FOR  RSCH  &  ENGR  R&AT  E  LS 
PENTAGON  ROOM  3D129 
WASH  DC  20301  3080 

NASA  SPACE  FLT  CTR 
ATMOSPHERIC  SCIENCES  DIV 
CODE  ED  41  1 
HUNTSVILLE  AL  35812 

NASA  MARSHAL  SPACE  FLT  CTR 
ATMOSPHERIC  SCIENCES  DIV 
E501 

ATTN  DR  FICHTL 
HUNTSVILLE  AL  35802 

NATIONAL  SECURITY  AGCY  W21 
ATTN  DR  LONGBOTHUM 
9800  SAVAGE  ROAD 
FT  GEORGE  G  MEADE  MD  20755  6000 

NAVAL  AIR  DEV  CTR 
CODE  5012 
ATTN  AL  SALIK 
WARMINISTER  PA  18974 

NAVAL  OCEAN  SYST  CTR 
CODE  54 

ATTN  DR  RICHTER 
SAN  DIEGO  CA  92152  5000 


NAV  RSRCH  LAB 

CODE  8150/SFA 

ATTN  J  BUISSON 

4555  OVERLOOK  DR  SW 

WASHINGTON  DC  20375  5354 

NAV  RSRCH  LAB 
ATTN  W  REID 
4333  OVERLOOK  DR  SW 
WASHINGTON  DC  20375 

NAV  SURFACE  WEAPONS  CTR 
CODE  G63 

DAHLGREN  VA  22448  5000 

NAV  SURFACE  WARFARE  CTR 
CODEK12 
ATTN  E  SWIFT 
17320  DAHLGREN  RD 
DAHLGREN  VA  22448  5100 

NIST 
MS  847.5 
ATTN  M  WEISS 
325  BROADWAY 
BOULDER  CO  80303 

NCAR  LIBRARY  SERIALS 
NATL  CTR  FOR  ATMOS  RSCH 
PO  BOX  3000 
BOULDER  CO  80307  3000 

ODCSOPS 

ATTN  BG  ROSE,  D  DSCHMIDT 
WASHINGTON  DC  20310  1001 

OIC  NAVSWC 

TECH  LIBRARY  CODE  E  232 
SILVER  SPRINGS  MD  20903  5000 
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ONR331 

ATTN  H  PILLOFF 
800  N  QUINCY  STRET 
ARLINGTON  VA  222 17 

OSD 

OUSD  (A&T)/ODDDR&E(R) 

ATTN  J  LUPO 
THE  PENTAGON 
WASHINGTON  DC  20301  7000 

OVERLOOK  SYSTEMS 
ATTN  D  BROWN 

1150  ACADEMY  PARK  LOOP  STE  1 14 
COLORADO  SPRINGS  CO  80910 

OVERLOOK  SYSTEMS 
ATTNTOCVIRK 

1150  ACADEMY  PARK  LOOP  STE  1 14 
COLORADO  SPRINGS  CO  80910 

PACIFIC  MISSILE  TEST  CTR 
GEOPHYSICS  DIV 
ATTN  CODE  3250 
POINT  MUGU  CA  93042  5000 

PAQ  COMMCTN 
ATTNQHUA 
607  SHETLAND  CT 
MILPTAS  CA  95035 

PHILLIPS  LABORATORY 
PLLYP 

ATTN  MR  CHISHOLM 
HANSCOM  AFB  MA  01731  5000 

PHILLIPS  LABORATORY 
PL  LYP  3 

HANSCOM  AFB  MA  01731  5000 
PL  WE 

KIRTLAND  AFB  NM  87 1 1 8  6008 


PRIN  DPTY  FOR  ACQUISITION  HDQTRS 

US  ARMY  MAIL  CMND 

AMCDCG  A 

ATTN  D  ADAMS 

5001  EISENHOWER  AVE 

ALEXANDRIA  VA  22333  0001 

PRIN  DPTY  FOR  TECHLGY  HDQTRS 
US  ARMY  MATE  CMND 
AMCDCG  T 
ATTN  M  FISETTE 
5001  EISENHOWER  AVE 
ALEXANDRIA  VA  22333  0001 

ROCKWELL  CACD 
ATTN  L  BURNS 
400  COLLINS  RDNE 
CEDAR  RAPIDS  lA  52398 

ROCKWELL  COLLINS 
ATTN  C  MASKO 
400  COLLINS  RDNE 
CEDAR  RAPIDS  lA  52498 

ROCKWELL  DA85 
ATTN  W  EMMER 
12214  LAKEWOOD  BLVD 
DOWNEY  CA  92104 

ROCKWELL  SPACE  OPS  CO 
ATTN  B  CARLSON 
442  DISCOVERER  AVE  STE  38 
FALSON  AFB  CO  80912  4438 

ROCKWELL  SPACE  OPS  CO 
AFMC  SSSG  DET  2/NOSO/ROCKWELL 
ATTN  R  SMETEK 
442  DISCOVERER  AVE  STE  38 
FALSON  AFB  CO  80912  4438 
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ROCKWELL  SPACE  SYSTEMS 
MAILCODE  841  DA49 
ATTN  D  MCMURRAY 
12214  LAKEWOOD  BLVD 
DOWNEY  CA  90241 

SCI  AND  TECHNOLOGY 
101  RESEARCH  DRIVE 
HAMPTON  VA  23666  1340 

SP  &  TERRESTRIAL  COMMCTN  DIV 
AMSEL  RD  ST  MC  M 
ATTN  H  SOICHER 
FT  MONMOUTH  NJ  07703  5203 

SPACE  ENVIRONMENT  LAB/NOAA 
R/E/SE 

ATTN  J  KUNCHES 
325  BROADWAY 
BOULDER  CO  80303 

SPECIAL  ASSIST  TO  THE  WING  CMNDR 
50SW/CCX 

ATTN  CAPT  H  BERNSTEIN 
300  O’MALLEY  AVE  STE  20 
FALCON  AFB  CO  80912  3020 

SRI 

M/S  BS378 
ATTN  M  MOEGLEIN 
333  RAVENSWOOD  AVE 
MENLO  PARK  CA  94025 

STANFORD  TELECOM 
ATTN  BF  SMITH 
1221  CROSSMAN  AVE 
SUNNYVALE  CA  94088 


STANFORD  UNIVERSITY 
HELP/GP  B 
ATTN  D  LAWRENCE 
STANFORD  CA  94305  4085 

STANFORD  UNIVERSITY 
HEPL/GP  B 
ATTN  T  WALTER 
STANFORD  CA  94305  4085 

TASC 

ATTN  T  BARTHOLOMEW 
1 190  WINTERSON  RD 
LINTHICUM  MD  21090 

TRIMBLE  NAV 
ATTN  L  KRUCZYNSKI 
585  N  MARY 
SUNNYVALE  CA  94086 

TRIMBLE  NAV 
ATTN  P  TURNEY 
585  N  MARY 
SUNNYVALE  CA  94086 

USAASA 

MOAS  AI 

ATTN  W  PARRON 

9325  GUNSTON  RD  STE  N3 19 

FT  BELVOIR  VA  22060  5582 

USAF  ROME  LAB  TECH 
CORRIDOR  W  STE  262  RL  SUL 
26  ELECTR  PKWY  BLD  106 
GRIFFISSAFB  NY  13441  4514 

USAF  SMC/CED 
DMA/JPO  ATTN  M  ISON 
2435  VELA  WAY  STE  1613 
LOS  ANGELES  AFB  CA  90245  5500 
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USAFETAC  DNE 
ATTN  MR  GLAUBER 
SCOTT  AFB  IL  62225  5008 

US  ARMY  CECRL 
CECRL  GP 
ATTN  DR  DETSCH 
HANOVER  NH  03755  1290 

USATRADOC 
ATCD  FA 

FT  MONROE  VA  2365 1  5 1 70 
USNO 

ATTN  B  BOLLWERK 
5580  PIEDRA  VISTA 
COLORADO  SPRINGS  CO  80908 

WSMR  TECH  LIBRARY  BR 
STEWS  IM  IT 
WSMR  NM  88001 


Record  Copy 


TOTAL 


